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October, 1922, Bulletin of the National Research Council. 
there is a table giving the mass absorption coefficients of a number of 
elements from Li to Pb for wave-lengths from .1A to 1.0A. For the most 
part these results are taken from the works of Hewlett,' Richtmyer,? 
and Duane* and are for homogeneous rays from a tungsten target 
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THE ABSORPTION COEFFICIENTS OF HOMOGENEOUS 
X-RAYS BETWEEN WAVE-LENGTHS 0.1 A AND 0.71 A 


By S. J. M. ALLEN 
ABSTRACT 


Using the general radiation from a tungsten target, dispersed by reflection 
from a crystal, the total mass absorption in fourteen elements of atomic number 
6 to 83 (C, Mg, Al, S, Fe, Ni, Zn, Ag, Sn, W, Pt, Au, Pb and Bi) was deter- 
mined for various wave-lengths from .10 A to.51 A by the ordinary ionization 
measurements. The Ka and K@ radiation from a Mo target was used to give 
values for wave-lengths .63 A and .71 A. Very good agreement was found 
with the results of Richtmyer and Duane wherever these are comparable. For 
the elements from C(6) to W(74) the results are represented within a few per 
cent by the formula p/p=1/p+o/p=Cr? *N?+40/p, where C=7.82(10)-* 
for the K absorption and 1.00(10)-* for the L absorption and the scattering 
coefficient ¢/p is assumed to increase with N from .144 for Al to .50 for W at 
.12 A, and also with the wave-length. The jump in the value of the fluorescent 
absorption coefficient 7/p from L to K absorption is the ratio of Cx to Cr or 7.8. 
The atomic absorption 7, =2.18(10)-*A2**N* (Bragg-Owens law). For the 
heavier elements Pt, Au, Pb and Bi the values of Cx get progressively lower to 
4.7(10)-* and those of Cy increasingly higher to 1.12(10)-*; the ratio or jump 
decreasing to 4.1 for Bi. 


INTRODUCTION 


N excellent summary of our knowledge of the absorption of x-rays 
when they pass through matter is given by A. H. Compton in the 


reflected from a calcite crystal. 
A careful examination of these results shows a considerable divergence 





of the mass absorption given by the different observers for the same 


1C, W. Hewlett, Phys. Rev. 17, 284 (1921) 
* Richtmyer, Phys. Rev. 18, 13, (1921) 
3 W. Duane, Proc. Nat. Acad., March 1922 


On page 32 
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substance and same wave-lengths. In the case of aluminum those of 
Hewlett are uniformly lower than those of Richtmyer. Those of Duane, 
over the short range investigated, agree well with those of Richtmyer. 
Only two results for the K absorption of an element of higher atomic 
number than that of silver are given, those of Hull and Rice‘ for lead, 
obtained eight years ago. It therefore seems desirable, considering the 
importance of the subject, to have a more complete set of absorption 
coefficients, comprising a greater range of wave-lengths and atomic 
numbers, especially the K absorption coefficients of the high atomic 
number elements. Experiments for this purpose have been carried out 
by the author during the last two years, and a summary of the results 
was given at the meeting of the American Physical Society at Cincinnati, 
December 1923. This paper gives the results and conclusions in a more 
detailed manner. 

If a narrow beam of homogeneous x-rays passes through a layer of 
any substance whose thickness is d and density p the loss of energy is 
given by, Ia4=IJoe“4 where Jy represents the energy on entrance, and 
Iq the energy on exit. 

The quantity yu/p=(log I)—log Ja)/pd, is known as the mass ab- 
sorption coefficient. From y/p can be calculated the atomic absorption 
coefficient u, by use of Avogadro’s constant and the atomic weight. 

The loss of energy represented by u/pis due to the energy absorbed to 
produce the various fluorescent radiations represented by r/p, and the 
energy lost by scattering etc. represented by o/p. 

Hence, u/p=r/p+a/p 
and, m ae Ma=Tat@a. 
It has génerally been accepted that these absorption coefficients can be 
represented empirically by the following expressions: 
t/p=Cr"N? 
Ta=C,r\" N¢ 
where \ is the wave-length in angstrém units; N the atomic number; 
Cand C, constants for each fluorescent series; 1, p, g, empirical exponents. 


Since »/p alone can be measured by experiment it is necessary in some 
manner to assume values of o/p, before +/p can be estimated. In the 
range of \ from .1 A to 1.0A, o/p in general is a quite large part of u/p, 
and for short wave-lengths and low values of N it is indeed the greater 
part, hence the uncertainty of o/p renders confirmation of the formulas 
difficult. The most usual values accepted for the exponents are, n=3, 
p=nearly 3, and q=4. 


‘ Hull and Rice, Phys. Rev. 8, 836 (1916) 
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When for any substance a wave-length is reached which is very close 
to its emission wave-length Ky, the absorption coefficient suffers a very 
abrupt change, becoming 8 times as great in the case of passing from 
LtoK. If the value of (7+/p)x for the K absorption just after this change 
takes place is known, and the value of (7p), for the same wave-length 
which would have been observed if this change had not taken place is 


also known, we can divide (7/p)x by (7/p), and obtain what we may call 


the “jump” J from L to K series. If the above formulas are sufficient 
to represent the values of (r/p)x and (r/p)z over the complete possible 
range of N, this jump should always be the same, about 8, and con- 
sequently Cx and C, should be constants. 


EXPERIMENTAL ARRANGEMENTS 


The modern method of measuring the absorption coefficients of 
homogeneous rays is well standardized, and it is not necessary to go into 
details here. The present work follows closely the methods used by 
Richtmyer and Duane. 

A fine beam of the general radiation from a Coolidge Universal Tung- 
sten tube was reflected from the (100) planes of a selected crystal mounted 
on a Bragg x-ray spectrometer. Both calcite and rock-salt crystals 
were used. 

The resulting narrow spectrum was intercepted by a fine slit in the 
ionization chamber which was filled with methyl iodide vapor. Readings 
were carried out as follows. The zero position of the spectrometer was 
carefully determined by setting on the Ka lines of Tungsten and also 
on the Ka and Kf lines of a Mo water-cooled tube. The tube was then 
excited to the desired voltage and the ionization chamber set to receive 
the definite wave-length to be tested. From the chamber angle 26, and 
the crystal spacing d, the value of \ was calculated from Bragg’s classical 
equation NA=2dsin 6. 

lonization readings were then taken both before and after an absorbing 
layer was placed in position between the target and the crystal. For 
any wave-length chosen, a complete absorption test to extinction was 
always made on some substance, such as aluminum or copper. By 
plotting the log J against the thickness of layer, a resulting straight line 
would ensure that all corrections had been made, and one would be 
sure that he was dealing with a pure homogeneous radiation. Then, 
keeping this wave fixed in position and the voltage constant, measure- 
ments were made on all the other elements. 

This method was used with all the wave-lengths observed from 

=.1A to \=.5A. The measurements at \=.63A and \=.71 A were 
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made with the Mo tube the a line being observed in the first, second, and 
third orders, and the 8 line in the first and second orders. 

Instrumental errors should not be greater on the average than 1 per 
cent. The chief difficulty, as has also been mentioned by Richtmyer, 
is to obtain very thin layers of the high atomic weight elements which 
are homogeneous. In all cases the mass of matter traversed was obtained 
by weighing a layer of known area. 

RESULTS 


The values of u/p obtained in these experiments are shown in Table 
I. In Fig. 1 are plotted, the results for the light elements C to S, in Fig. 
2 those for the heavier elements Fe to Ag. 


TABLE I 
Mass absorption coefficients p/p 
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These results cover a range of X N which is about the same as 
investigated by Richtmyer, Hewlett, and Duane. There is extremely 
good agreement between the results of the author and those of Richt- 
myer, and also those of Duane in the case of Al, Cu, and Ag. The 
results for C agree well with those of Hewlett at short wave-lengths, 
but are generally higher than his at longer wave-lengths. All elements 
show a steady but slower decrease as \ approaches .1A. There is no 
indication as yet of any abrupt change. In Figs. 3 and 4 are plotted 
those portions of the results for the high atomic number elements W, 
Pt, Au, Pb, and Bi between .10 A and .24 A, which include the critical 
absorption wave-lengths Ky. It is to be noticed that the observations 
are carried right through the abrupt changes of u/p. 
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These results for high weight elements show very clearly a distinct 
difference from similar results for lower weight elements. For example, 


if we compare the results for W and Bi, we see that in the L absorption 
60 


nA 
10 “30 
Fig. 1. n/p as a function of wave-length for C, Mg, Al and S. 


‘10 0 


Fig. 2. u/p as a function of wave-length for Fe, Cu, Zn and Ag. 


the graphs for Bi lie above that for W in accordance with the order of 
the atomic numbers, but in the K absorption this is reversed. The 
same phenomenon is true for all the heavy elements. In short, the sudden 
rise, or “jump” from’L to K absorption is much less than would be 
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predicted by extrapolation of the formulas given above, and moreover 
this jump becomes progressively less as N increases. 

In Table II are shown these “jumps” J calculated from the graphs 
and corrected for scattering from assumed values of o/p. 

The values of (7/p)x are the highest K values experimentally obtained 
and } is the wave-length at which they occurred. The values of (r/p)z 
are those at the same wave-length of the L absorption which would have 
been obtained if the “jump” had not occurred. With the high weight 
atoms are included also Ag, Mo (from Richtmyer) and Cu (from Barkla). 


TABLE II . 
Relations between the K and L absorption coefficients at the absorption limit Ky 








Element ry (u/p)e . a/p (r/p) x (u/p)y (r/p)y P 


1.46A 440 — 440 55 55 
.60 81 , 80. 11. 10.6 
.45 53.5 , 53. 7. 7.0 
.159 9.00 : . ; 1.40 
. 138 7.35 . 6: : , 1.27 
. 136 6.65 : : ; 1.21 
.125 $.25 i . ‘ 1.01 
.120 4.40 ; ; ; .90 
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It is readily seen that, while this ‘‘jump”’ in the absorption coefficients 
in passing from the L to the K levels is approximately 8 times for light 
elements, it becomes progressively less as N increases, showing a quite 
rapid decline from N=74 to N=83. 

The values of J here given depend upon the values of ¢/p assumed and, 
since we have no very exact knowledge of o/p for high weight atoms, 
they can be considered at this stage as only approximate. There seems, 
however, no doubt that they progressively decrease with the increase 
of N. 

Whether the absorption coefficients change so abruptly as to constitute 
a discontinuity is open to question. The results here indicated with high 
atomic number elements would seem to show that the abrupt rise slants 
backward towards shorter wave-lengths and that actual measurements 
can be made in this changing region, even with slit widths of .1 to.2 mm. 

At the same meeting mentioned above, Richtmyer read a paper on 
the relations of the absorption coefficients in the immediate neighbor- 
hood of the K critical absorption, which showed results almost identical 
with those of the author, confirming the idea that the value of J decreases 
steadily with the increase of N. 

A very careful analysis of all the values of u/p in this work was made 
to test the accuracy of the empirical formulas discussed above. This 
was done by actual calculation by trial and error instead of the method 
usually employed of plotting »/p against N*. It was found that, if 
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suitable values of o/p were chosen, over a wide range of \ and N up to 
47, the fluorescent absorption coefficient +/p could be expressed by a 
formula of the form, 
t/p=CrN?, 
and the atomic coefficient by 
Ta=C_rA"N’. 
a/p was taken for Al at \=.12A as .144 and at A=1.0A as .20, with 
increasing values as N increased. ; 
It was found that on the average n= 2.92 (instead of the usual value 
3), p=2.92, and q=4 best fitted the results. 
TABLE III 


Empirical data with reference to r/p= Cr?" N?-@ and ta = Cad? N 4 


es CK CL (Ca)K (Ca)L 
Element \ A=.12A (<10-3) (x10-3) (10°) (x10-%) 
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Extrapolation of these formulas to high atomic weight elements for 
the K series gave predicted values much higher than those actually 
measured. The accord was better in the case of the L absorption. The 
complete data are shown in Table III. The results for Mo are taken 
from Richtmyer and those for Sn from Bragg and Pierce. The universal 
constant (C,)x is=2.18  10-% for N=2.92 (K absorption) and (C,); 
=.29 X 10-* (L absorption). Cg =7.82 X 107°, and Cp=1.0 X 10°. 

It would seem that these formulas do not hold throughout the whole 
range of \ and N, and certainly not near a critical absorption region. 

Very notable attempts to deduce a theoretical formula for r. have 
been made recently by De Broglie’ and by H. A. Kramers® on the quantum 
theory. Both of these formulas indicate a decrease to be expected in 
J as N increases, which would be in accord with the results of this paper. 

UNIVERSITY OF CINCINNATI, : 

January 25, 1924. 


5 L. de Broglie, Jour. de Phys. et Rad. 1922 
*H. A. Kramers, Phil. Mag. Nov. 1923 





K AND L ABSORPTION SPECTRA OF TUNGSTEN 


THE K AND L ABSORPTION AND EMISSION 
SPECTRA OF TUNGSTEN 


By C. B. Crorutt 


ABSTRACT 


Both the absorption and emission spectra were photographed simul- 
taneously on the same plate, using a special float-operated mechanism to turn 
the crystal: very slowly so as to obtain the faint lines between ®: and 8;. The 
absorption \wave-lengths were found to be: LA,, 1.2122A; LA,, 1.0716A; LAs, 
1.0217A; KA, .17802A. Table I gives the emission wave-lengths of 28 L emission 
lines including the four new lines Bi; (1.2432), Big (1.2166), ye (1.0748) and 
n (1.0699), and of the Ks line (.18525) which is the line found by DeBroglie. 
Comparison with the Bohr theory. The agreement of the emission frequencies 
computed from the energy levels with the measured wave-lengths is better than 
1/10 per cent, except for 6:5 and Bis. The lines Bio, 89,712, and yu, however, cor- 
respond to transitions between energy levels which are not predicted by the 
selection principle, and lines 8s and n are not in agreement with the accepted 
energy levels. While Ln may be due to an impurity, #s is probably a tungsten 
line. 


HIS work is an investigation of the absorption and emission spectra 

of tungsten under the same experimental conditions. Both spectra 
were obtained simultaneously on the same photographic plate in order 
to make certain the relative positions of the absorption limits and the 
emission lines. 


EXPERIMENTAL METHOD 


A medium focus tungsten target Coolidge tube was used as a source 
of x-rays. The x-rays were analyzed by means of a rock salt crystal and 
the spectra were recorded on photographic plates. The apparatus is 
essentially that used by Dershem! with the substitution of an accurately 
constructed crystal table, a narrow slit of 0.005 cm width, and a new 
rotator for the crystal. 

The crystal table is shown in Fig. 1 and consists of a spindle fitted 
in bearings as shown. This spindle is removable and additional spindles 
make it possible to change crystals without remounting each time. 

The crystal was rotated through a small angle by means of levers 
attached to a float in a tank of water, the level of which was raised by 
means of a small stream of water and lowered by the use of an inter- 
mittent syphon. The advantage of this method is that the crystal can 
be set for a particular line and a prolonged exposure taken with the as- 


1 Dershem, Phys. Rev. 11, 461 (1918) 
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surance that the crystal will not rotate too far. The rotation can be made 
as slow as desired and the angle kept as small as a quarter of a degree. 
The fine slit greatly prolongs the time of exposure but enables greater 
resolving power to be obtained. The fine slit and the rotator made it 
possible to obtain the group of faint lines between #2 and £5. 

While working with the L series it was necessary to keep the voltage 
so low that no appreciable second order 
continuous radiation would affect the 

— plate at the same point. 

The work on the K series was done in 
the third order using 110,000 volts, 
effective value. The first and second 
order radiations were removed by ab- 
sorbing screens. 

The absorption bands were obtained 
by placing a screen of tungsten oxide in 
the path of the beam before it entered 
the spectrometer. This in no way inter- 
feres with the production of the emis- 
sion spectra but does interfere with the 
Bearing analysis as emission lines are close to 
the absorption bands. 

The method of measuring the plates 
is that developed by Dershem! in which 
a correction is made for the width of 
the slit. Measurements of faint emission lines and absorption bands 
were always made relative to strong emission lines. 














Circular Scale 











Ramune 
Fig. 1. Crystal table. 


CORRECTIONS 


With the method employed there are two chief sources of error. 

(a) If the plate is not perpendicular to the line of collimation it can 
be shown from geometrical considerations that 6 is given by the follow- 
ing formula 





D tan 20=dy/(1+m?) + /d?(1+m?)+Dm?, 


where 6 is the glancing angle, d the distance from the crystal to the plate, 
D the distance between the lines as photographed on both sides of the 
line of collimation, and m the tangent of the angular displacement of 
the normal to the plate from the line of collimation. Since this cor- 
rection is small, m may be determined with sufficient accuracy by means 
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of two known lines. In a manner similar to the above it can be shown 
that 





2(1—p)m = (1+ p)(mi+me2) + V(1+p)2(m +m) —4(1 — p)*myme 
where p is the ratio of the distances between the two lines on opposite 
sides of the line of collimation, m; and me the value of tan 6 for the two 
known lines. 

(b) If the face of the crystal is not in the axis of rotation an amount 
equal to 6/sin @ must be added to or subtracted from the distance 
from the axis of rotation to the plate in order to obtain the distance from 
the crystal to the plate, where 6 is the displacement of the face of the 
crystal from the axis of rotation. Since for small values of 6, sin@ is 
small this correction is important. 


TABLE I 
Wave-lengths in Angstroms 


L SERIES 
. 2432 Bio 
.2421 By 
. 2364 5 


1.4844 
1.4733 
1.418 


. 2099 
.2027 
. 1299 


1.2875 
1.2793 
1.2602 


.0786 
.0748 
.0723 


.2166 8 
. 2133 7Yi2 
.2122 6 


K SERIES 
Comparison lines Measured 
. 21352 BB .18436 Bs .18525 A .17802 
. 20885 y  .17940 


1 1 

1 1 

1 1 
1.2988 1.2217 vy 1.0964 

1 1 

1 1 

1 1 


Using Dershem’s method for locating the lines, the penetration into 
the crystal does not introduce an error for the method depends on locat- 
ing the edge of the line as reflected from the face of the crystal. 


EXPERIMENTAL RESULTS 

Table I shows the experimental results on the emission and absorption 
in the L and K series. The absorption limits are indicated by the letter A 
with appropriate subscript. In the L series the following lines were meas- 
ured carefully and the results taken as standards for the measurements 
of all other lines: a1, 81, 82, ¥1, ys, and ys. The comparison lines in the 
K spectrum are indicated in Table I. Because of the presence of the 
reference lines the results obtained on two or more plates varied only 
in the last significant figure shown. The line | does not appear because 


of the absorption in the glass walls of the tube. The lines Bis, Bis, Y12, 
and n are new lines for the L series of tungsten. 
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DISCUSSION OF RESULTS 


Fig. 2 presents the results in the well known form of energy levels as 
recently published by Bohr and Coster.? Table II gives the corresponding 
frequencies of the energy levels. The K and L frequencies were obtained 
directly from the measurements made on the wave-lengths. The other 


0 
0; 5(3,3) 
O2 5(3,2) 
0; 5(2,2) 
04 5(2,1) 
5(1,1) 


M 
3(3,3) 
3(3,2) 
3(2,2) 
3(2,1) 
3(1,1) 


2 fs 
a 
- . 


BALY 42 A4%4%19,0, 8,4 8A58,4%), 


K 
1(1,1) 


' 

' 
a, 

' 

' 


aa8A 7 
LET) 
Fig. 2. Energy level diagram for tungsten showing the K and L series. Solid 


lines are transitions predicted by the selection principle. Broken lines are other 
transitions. 


wal 





frequencies were calculated from the K and L absorption frequencies and 
the frequencies of the emission lines. The solid lines between the energy 
levels of Fig. 2 represent transitions predicted by the selection principle.’ 


? Bohr and Caster, Zeits. f. Phys. 12, No. 6 (1923) 
3 Bohr, Phil. Mag. June 1922 
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The broken lines represent other transitions found. The notation is that 
of the line produced by the transition. If there is no symbol the line has 
not been found. The symbols at the left of the horizontal lines are the 
usual notations of the energy levels; those at the right are the notations 
of Bohr and Coster [the symbol 4(2,1) corresponds to the energy level 
having the total quantum number 4, the “grund’’ quantum number 2 
and the azimuth quantum number 1]. 

The selection principle states that the azimuth quantum number must 
change by +1 or 0 and that the “grund”’ quantum number must change 
by +1.? 

TABLE II 
Absorption frequencies X 10°" 
K and L measured; M, N and O calculated 


K, 16.844 M, .4401 .0579 .0022 


Li 2.4738 M; .4548 .0629 .0053 

L, 2.7982 M; .5554 . 1058 .0140 

L; 2.9349 My, .6261 . 1220 .017 
Ms .6814 . 1448 .0189 


A comparison of Table I with Fig. 2 shows that all but the 6s and n 
lines are represented on the diagram. The n line is a new line and may 
be due to some impurity, although a comparison with tables failed to 
reveal its identity. The 8s line has been found by Overn‘ and Siegbahn’* 
and should be accepted as a line of tungsten. This line is hard to measure 
being very indistinct. Coster® has shown that it appears as B,; and Bj. 
in elements of atomic numbers 47 to 51. Of the four lines of the L series 
Bo, Bio, Yi1, and 712 not given by the selection principle, 712 is a new line. 
A careful investigation failed to reveal the Ka; line of which Duane and 
Shimizu’ found some indication. The results indicate that if it does exist 
it does not differ in wave-length from the az line by more than .04 per cent, 
or else its intensity is so low that it is concealed by the general radiation. 
It is not predicted by the selection principle. The two other new lines 
are given by the selection principle. 

Bis is so close to B, that it could not be completely resolved, but be- 
cause of the difference in intensity the edges of both lines could be located 
and measured. The other new lines are very faint and close to lines 
already known. Kf; is the line previously found by DeBroglie.* The 
wave-lengths obtained for the absorption bands are shorter than those 


* Overn, Phys. Rev. 14, 137 (1919) 

5 Siegbahn, Phil. Mag. Nov. 1919 

® Coster, Phil. Mag. Sept. (1922) 

7 Duane and Shimizu, Phys. Rev. 14, 67 (1919) 
® De Broglie, Comptes Rendus, 1920 
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obtained by Duane, which are usually accepted as standard. These 
changes in wave-length are sufficient to alter the relative positions of 
the LA; and the LA, with respect to 8; and y¢, as is shown in Table III 
and of course visually on the plates themselves. 

It is seen from Fig. 2 that all lines found as transitions between energy 
levels but not predicted by the selection principle involve (1,1) energy 
levels. This may be a coincidence but it may have significance. 

No special attempt was made to obtain the relative intensities of the 
lines. The 61; is of medium intensity while the other new lines are very 
faint. 

The fact that the lines 85, Bis, and y12 fit into the energy level diagram 
makes it seem reasonable that they are tungsten lines. 

Since the absorption and emission spectra appear on the same plates 
it is impossible to say much regarding the structure of the abscrption 
bands, some of the emission lines being very close to the absorption edges. 
However the edges were sharp and easily: measured. 


TABLE III 


The shift in absorption wave-lengths 


Ai 1.2140* A: 1.0730* ys 1.0266 
B, 1.2133 ys 1.0723 A; 1.024* 
A; 1.2122 A; 1.0716 A; 1.0217 


The differences between the frequencies of the emission lines as com- 
puted (1) from the energy levels and (2) from the measured results are 


about .13 per cent for 8:; and Bi. and less than .10 per cent for the re- 
mainder with an average difference of .03 per cent. 


CONCLUSIONS 


The selection principle reduces the possible number of lines from fifty- 
one to twenty-one for the L series and from twenty to eight for the K 
series. Two of the new lines found agree with the selection principle, thus 
increasing the total number in agreement from eighteen to twenty for 
the L series, leaving one more line to be found. All lines in the K series 
are in agreement with the princip!e. 

In addition to the lines given by the selection principle four were found 
in the L series which are not in agreement but which quantitatively 
correspond to transitions between energy levels. The Ka; line, found by 
Duane but not observed by the present photographic method, is not in 
agreement with the selection principle and is yet somewhat in doubt. 

As an indication that the present energy levels are not sufficient there 
appear to be two lines in the L series which are not in agreement. One 


*Values given by Duane and Shimizu’ 
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of these is a new line and may be due to an impurity; the other is the 
Bs line. 

It is seen from the above that the energy levels as drawn do not satisfy 
all the measurements of what are undoubtedly tungsten lines, and that 
Bohr’s selection principle would eliminate four lines which have been 
found. These disagreements emphasize the need for further experimental 
and theoretical work. 

I wish to acknowledge indebtedness to the members of the staff of 
the physical laboratory at the University of Iowa for their interest in 
the work and especially to Professor G. W. Stewart who suggested the 
problem. 


UNIVERSITY OF IOWA. 
February 27, 1924. 








| 
| 
| 
| 








J. J. NOLAN 


THE CONSTITUTION OF GASEOUS IONS 
By J. J. NoLAN 


ABSTRACT 


Mobilities of various groups of ions in dry and moist air.—Using the 
Franck modification of the Rutherford alternating current method, with the 
plate and gauze accurately parallel and 2.5 cm apart, current-voltage curves 
were obtained as the 52 cycle voltage was increased from 200 to 500 volts. 
Previous experiments have shown the presence of four predominating types 
of positive ions with mobilities of 2.04, 1.73, 1.52 and 1.34. The curves now 
reproduced show periodic inflections about every 10 volts which indicate 
intervening groups making some twenty groups of ions with mobilities varying 

_from 2.24 to 1.34. These values agree closely with those computed according 
to J. J. Thomson's theory for clusters of from 15 to 36 water-molecules. In dry 
air the groups with the lower mobilities are less prominent. Other experi- 
mental evidence is cited in favor of the cluster theory as opposed to the small 
ion theory of ions. Loeb’s criticism of the author’s work is discussed and it is 
shown that there is no real divergence between these results and those of Loeb. 


HE constitution of the ions in gases at normal pressures is still a 

matter of controversy. In various papers’** the present writer has 
brought forward evidence which, if accepted, would be decisive in favor 
of the older cluster theory as against the more recent view which regards 
the gaseous ions as simply charged molecules. It was claimed that these 
experiments demonstrate that the normal ionization consists of a mixture 
of distinct groups having mobilities not widely different from one another, 
and that in addition there are present in small numbers groups of ions 
of much higher mobility. These conclusions have been contested by 
Loeb‘ who holds that his experiments demonstrate the existence of only 
one type of ion. He is of opinion that results which suggest a complex 
ionization are spurious and are due to “failure to observe some of the 
fundamental conditions’’ of successful experiment. The object of the 
present paper is, firstly to describe fully the apparatus and the experi- 
mental method used by the author so as to permit of criticism on that 
head being definite; secondly to show that the results of Loeb when 
properly interpreted, not merely do not controvert but actually support 


1 J. J. Nolan, Proc. Roy. Irish Acad. 35, 38 (1920) 

?J. J. Nolan, Proc. Roy. Irish Acad. 36, 74 (1923) 

*J. J. Nolan and J. T. Harris, Proc. Roy. Irish Acad. 36, 31 (1922) 
*L. B. Loeb, Journal Franklin Inst. 196, 537 (1923) 
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the contentions of the present writer; and finally to bring forward further 
evidence as to the structure of ions in air. 


APPARATUS AND EXPERIMENTAL METHODS 


The method of measuring ionic mobilities employed in the work under 
discussion is the well-known modification of Rytherford’s method due 
to Franck. lons are produced on one side of a metallic gauze and those 
of one sign are driven by a steady field through the gauze. Here they 
come under the influence of an alternating field of sine form between the 
gauze and a parallel plate. If the potential difference is such that its 
value read on an electrostatic voltmeter is greater than rnd?//2u where 
n is the frequency, d the distance between the gauze and plate and u 
the mobility of the ions, some of the ions coming from the gauze will be 
captured by the plate before the reversal of the field. For fields less than 
this value the plate will receive no charge. If the current to the plate is 
plotted against the voltage and if there is only one kind of ion present, 
a smooth curve should be obtained. This curve should intersect the volt- 
age axis sharply at the critical voltage, and from the value of this the 
mobility of the ion can be readily determined. If there is more than one 
kind of ion present, the slope of the curve will increase sharply at the 
critical voltage peculiar to each type of ion. 

As the apparatus employed has varied somewhat from time to time, 
it is proposed to describe the particular form with which the results given 
later in this paper were obtained. The modifications adopted in other 
work will then be indicated. The essentials are shown in Fig. 1. 

The plate which receives the ions was mounted as follows. A thick 


brass plate, 13 cm in diameter was securely screwed on to an ebonite 


plate of the same size and 1.3 cm thick. The face of the plate was then 
turned flat in the lathe and a circular cut less than 1 mm in width was 
made separating an inner circular plate 8 cm in diameter from a guard- 
ring, which was thus a little less than 2.5 cm wide. The plate system was 
supported by a stout brass rod passing through an ebonite plug in the 
iron plate which formed the base of the apparatus. This rod served as a 
connection between the insulated plate and the electrometer. The guard- 
ring was connected to the iron plate and thus to earth. The gauze, 8 cm 
in diameter, was fitted to a brass ring of the same dimensions as the 
guard-ring. The gauze was made of brass wire .35mm in diameter and 
had approximately 50 meshes per cm?. Great care was taken in soldering 
the gauze to the brass ring to make it quite taut so that no sagging could 
occur. It was necessary to insure that the lower surface of the gauze was 
flush with the ring. A flat circular plate was prepared slightly less than 
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8 cm in diameter and of a thickness equal to the difference between the 
thickness of the ring and of the gauze. This plate was placed against the 
gauze on the upper side and the whole clamped firmly in a vise over night. 
In consequence of these precautions, the under surface of the gauze was 
quite level, and its distance from the plate was given very accurately 
by the height of the three ebonite pillars supporting its rim. The pillars 
were carefully turned in the lathe. Their diameter was 1.1 cm and their 
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Fig. 1. Apparatus for determining? mobilities by the Rutherford-Franck methcd. 








mean height 2.494 cm, the three differing only slightly in the last figure. 
The pillars stood at the outer edge of the guard-ring, and a minute trace 
of seccotine was used to fasten them in position. The rim of the gauze 
rested on them without any adhesive. 

The source of ionization was a deposit of polonium on the central 
portion of a small disk of bismuth which rested on the centre of the gauze. 
At a distance of 5 cm above the gauze was a parallel plate and between 
this and the gauze a steady potential difference was maintained by means 
of dry cells. In most of the work to be described subsequently, the value 
of this potential difference was 50 volts. The plate was fitted with a small 
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rim and was used to support porcelain dishes containing phosphorus 
pentoxide when drying was desired. The plate was supported by a rod 
passing out through an ebonite plug which served as a stopper in the 
glass bell-jar which surrounded the apparatus. A wire making connection 
with the gauze was brought out through the same plug. Part of the wire 
inside the vessel was surrounded by a thin glass tube. The bell-jar was 
16.5 cm in diameter internally. It was ground on to the iron base-plate, 


the contact being made air-tight by the use of a little vaseline. The 
ebonite plug in the neck of the bell-jar was sunk about } cm and the space 


flooded with molten paraffin wax. 

The alternating potential was derived from the secondary of a trans- 
former, one terminal of the secondary being earthed and the other con- 
nected to the gauze through a 2-volt cell. The cell was so connected that 
the upward field on the ions was slightly greater than the downward, 
so that ions which had just escaped capture at the plate were discharged 
by the gauze, or at least compelled to traverse the full upward distance 
before the next alternation of the field. The primary of the transformer 
was connected to the Dublin supply, the frequency of which is about 52 
The voltage of the secondary, which could be readily varied by varying 
the resistance in the primary, was read off directly on an Ayrton and 
Mather electrostatic voltmeter. Two of these were used; one having the 
range 80-250 volts could be read accurately to } volt, the other having 
the range 200-600 volts could be read accurately to 1 volt. The frequency 
was measured at the beginning and end of each set of observations and 
was remarkably steady. The electrometer used was of the Dolezalek 
type with quartz insulation and a phosphor-bronze suspension. Its 
sensitivity was about 1000 mm/volt and the capacity of the electrometer 
and apparatus was 165 cm. 

VARIATIONS IN APPARATUS 

Gauze. A number of sets of observations were made with a perforated 
brass plate substituted for the gauze. The plate was quite fat, 1 mm in 
thickness and had 12 circular holes per cm’, each 2.35 mm in diameter. 
The results did not differ in general character from those obtained with 
the gauze except that the currents to the lower plate were considerably 
smaller. It is evident that a smaller number of ions found their way 
through the perforated plate. 

Vessel. In about half of the observations a metal enclosure was used 
instead of the glass bell-jar: The dimensions were almost exactly the 
same except that the top was flat. In this case the connection to the gauze 
was brought out through the iron base-plate. As will be seen later, the 
essential character of the results was not affected by these modifications. 
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Apparatus for hydrogen. The apparatus used for observations on 
hydrogen,? which was also used for some work on air, differed from that 
described only in dimensions. The electrometer-plate was 10 cm in 
diameter and the guard-ring 15 cm. The height of the pillars was exactly 
4cm. The apparatus was enclosed in a brass vessel 20 cm in internal 
diameter. 


CONDITIONS FOR ACCURATE OBSERVATIONS 


Constancy of ionization. This is secured by the use of polonium. From 
some points of view, ions produced photo-electrically are preferable, since 
any uncertainty as to conditions at the grid is avoided. On the other hand, 
one is limited to negative ions, and, a more serious matter, steadiness of 
ionization is difficult to attain. 

Constancy of alternating field in amplitude and frequency. This is most 
likely to be attained from a large a.c. system such as that of the Dublin 
public supply. Except at certain hours of the day, the voltage as recorded 
by the voltmeters was remarkably steady. The practice was adopted of 
working during intervals in which previous experience had shown that 
the voltage was undisturbed. As far as may be judged from the volt- 
meter readings, the voltage on favorable occasions was steady to at least 
3 per cent. The curves obtained bear out this view. The maximum change 
of frequency noted from day to day was of the order of 1 per cent; the 
maximum variation between the beginning and end of an experiment 
was about 2 parts in 500. 

Parallelism of gauze and plate. The distance between the gauze and the 
plate enters into the formula for the determination of mobility in the 
second power. It is obviously of great importance that exact parallelism 
should be secured. The apparatus used by the author has been specially 
designed in order to make sure of this point. In the form of apparatus 
used by Loeb, strict parallelism is very difficult to obtain. Loeb estimates 
an uncertainty in his mobility values of +3 per cent on this account. 

Uniformity of field. The guard-ring should effectively eliminate any 
distortion of the field due to the walls of the vessel. The fact that the 
substitution of metal for glass made no difference in the results shows 


that the protection is effective. It readily appears from calculation that 
the density of ionization used will not appreciably disturb the field. There 
remains the possibility which has been suggested® of penetration of fields 
through the gauze. If the lines of force from the upper (accelerating field) 
penetrate through the gauze in any numbers, one would except to find 
a current to the lower plate when the gauze is at earth. An effect of this 


5 Miss Zimmerschied, Phys. Rev. 21, 721 (1923) 
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character has been noted previously*® and a different explanation sug- 
gested. This effect appeared only when the intensity of ionization was 
high and when the gas was very dry. It has not been observed in the 
course of the work dealt with in this paper. 

Induction effects due to the alternating field. Various devices are em- 
ployed by different workers in order to eliminate this cause of disturbance 
to the electrometer system. The simplest method is to disconnect the 
electrometer from earth before closing the primary circuit of the alternat- 
ing system. The inductive effect is then almost negligible and the elec- 
trometer needle soon takes up its uniform rate of motion. 

Insulation. It is essential that the insulation of the electrometer sys- 
tem should be good. The test of insulation applied was that the elec- 
trometer and plate should not leak more than 5 divisions per minute 
when charged to .14 volt, no condenser being in parallel. 


RESULTS OF EXPERIMENTS 


In Fig. 2 are plotted some of the current-voltage curves obtained with 
negative ions. Curve A is the result of some observations made on air 
saturated with water-vapor. The lower part of the curve was obtained 
with a condenser of capacity .001 mf in parallel with the electrometer, 


the upper part with a capacity of .003 mf. The readings with the higher 
capacity were multiplied by the appropriate factor to make them cor- 
respond to the readings with the lower capacity. The two sets of observa- 
tions overlap between 470 and 490 volts. It is not possible with any pre- 
tence to accuracy to draw through these points a smooth curve cutting 
the axis sharply. A nearer approximation represents it as a number of 
straight lines intersecting at 520, 474 and 435 volts. The constant 
(mobility X critical voltage) was 710. The curve therefore represents ions 
of mobilities 1.36, 1.50 and 1.63, corresponding to three of the ions to 
which the author attributes the bulk of the ionization in moist air, The 
curve shows the presence of another ion, but the intercept with the volt- 
age axis is not sharp enough to allow of its mobility being calculated. 
Now it is not contended that the curve given is a simple composite of 
four curves, one for each kind of ion. It will be seen later that close ex- 
amination would show it to be more complex. But it does show that the 
bulk of the ions in moist air tend, as has been claimed in previous papers, 
to fall into four principal groups of which three are here indicated. 

A remaikable confirmation of this result has been found in the cloud 
experiments of C. T. R. Wilson.’ In one photograph, the separation in 
an electric field of the four kinds of ion produced in the track of an a- 


6C. T. R. Wilson, Proc. Camb. Phil. Soc. 21, 405 (1923) 





22 J. J. NOLAN 


particle is shown by four parallel cloud-tracks. It is difficult to conceive 
of any possible explanation of this photograph other than the four-fold 
character of the ionization. 
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Fig. 2. Current-voltage curves 
A, saturated air; B, air partially dried; C and D, dry air; E, Loeb’s curve. 





In the work described in the present paper, the interest was not so 
much in the ionization in moist air as in the changes produced by drying 
and the accurate estimation of the mobilities of the faster groups of ions. 
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Porcelain dishes containing phosphorus pentoxide were inserted and the 
vessel was sealed up. Two days after the drying began Curve B was 
obtained. The upper part of this curve was taken with a capacity of 
.008 mf in parallel with the electrometer, the lower parts with .003 and 
.001 mf. The overlapping parts agree very well except for some observa- 
tions between 390 and 400 volts. The ordinates of this curve (and of 
Curves C and D) should be multiplied by 2.67 for comparison with 
Curve A. If what Loeb calls the ‘asymptotic foot” at the bottom of the 
curve be disregarded, a straight line may be drawn through the plotted 
points. This line intersects the voltage axis at 352, and it might be said 
that the observations therefore indicate the existence of a single ion of 
mobility 710/352 =2.02. But in the opinion of the writer, this summary 
treatment does not reveal the true conditions. In the first place the 
asymptotic feet in these curves cannot be so lightly disregarded; and in 
the second place, the straight-line part of the curve will, if attentively 
examined, show a remarkable peculiarity. It has a gapped appearance 
due to the fact that while the observed voltages are fairly uniformly 
distributed, the current values show sudden jumps at regular intervals. 
In the spaces between the gaps, the current values, while never very far 
from the straight line, exhibit a tendency to lie first above and then 
below it. In fact it would seem that the curve is not a straight line but a 
series of short curves, each short curve corresponding to a step of about 
10 volts. It might be alleged that this appearance is quite illusory; that 
the regularly occurring curvature is due to some experimental error— 
some error in reading the voltmeter for example. Now it will be shown 
later when the curve is examined closely and not in the rather crude 
fashion just given, that this complex structure can be quite definitely 
and accurately established. On the other hand, Curve C which was ob- 
tained with very dry air, is plotted to show that at a higher degree of 
drying the upper part of the curve is quite smooth. If the form of Curve 
B is due to experimental irregularities of any kind, they should equally 
well affect Curve C, for the only difference is that one curve was plotted 
40 days later than the other. 

In order to compare the results with those obtained by Loeb, Curve D 
has been plotted. This curve was obtained 23 days after drying had begun, 
that is 21 days after Curve B. All observations were taken with a con- 
denser of .008 mf in parallel with the electrometer. Curve E has been 
plotted from Loeb’s observations taken from Fig. 3 of his paper. The 
points have been taken from an enlarged photograph of Loeb’s curve 
and it is believed that they are accurately reproduced. While in general 
outline there is not much difference between the author’s curve and that 
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of Loeb, it will be noticed that the latter is not of as regular a character. 
Loeb draws a smooth line through his points, but the author has ventured 
as the diagram shows, to follow the points rather more closely. Possibly 
Loeb’s observations are inaccurate owing to want of steadiness in the 
ionization and want of parallelism is his apparatus, but as he states that 
each point was obtained from the average of two concordant electrometer 
readings, it may be assumed that inaccuracy on the first head at least 
has been eliminated. This being so, the author cannot find any reason 
for refusing to give to Loeb’s experimental points their full value, es- 
pecially as they can be made to lie accurately on a number of smooth 
curves instead of very inaccurately on one. There can be no doubt 
whatever that Loeb’s curve shows a bend at 115 volts. It also shows 
but not so strikingly, bends at about 83.5 volts and 78 volts and there is 
a slight indication of a sag in the curve at about 103 volts. Loeb’s mobility 
constant was 169 approximately. Therefore his curve shows quite dis- 
tinctly the ion of mobility 169/115 = 1.47, also ions of mobility 2.02 and 
2.17. It shows indicationg of a possible ion of mobility about 169/103 
= 1.64. It is hardly necessary to point out how closely these results agree, 
not merely in general character, but in the actual values of the mobilities 
with the results published by the present writer. The curve for dry 
air (D) is much smoother than that of Loeb, but seems to have the 
“asymptotic foot”? better developed. It is, in fact, what the author is 
accustomed to regard as the characteristic curve for a dry gas, while 
Loeb’s curve is intermediate in character between that for a dry gas and 
that for a saturated gas (A). In other words, the drying in Loeb’s case 
seems to have been imperfect. The question of drying is not a simple 
one, but the work of Baker’ and others seems to show that the gas in a 
vessel cannot be dried by simply filling it through drying materials. 
Prolonged contact of the gas with the drying agent is necessary. 

Loeb’s other curves (Figs. 1 and 2 of his paper) are not so well suited 
for comparison, principally because they show too small a number of 
points over the range under discussion. Inside that range, however, they 
certainly show irregularities of such a character as to suggest that if 
more points were determined the complex character of the curves would 
be disclosed. 

The Curves A to D have been plotted only to show that there is no 
divergence between the author’s results and those of Loeb so far as the 


latter has gone, and that it is not necessary to seek for explanations 
of differences which have not been proved to exist. Curves of this char- 


? Baker, Jour. Chem. Soc. (London) 81, 400 (1902); 101, 2339, (1912); 121, 568 (1922) 
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acter have been regarded by the author as only in the nature of prelimi- 
nary surveys of the ground. The real investigation consists in the close 
scrutiny of the current-voltage curve stage by stage. 
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300 3/0 520 S30 340 J50O 
Fig. 3. Current-voltage curves for 300-350 volts, at different stages of drying. 





ACCURATE INVESTIGATION OF THE STRUCTURE 
OF CURRENT-VOLTAGE CURVES 


In Fig. 3 a number of curves are plotted showing the examination of 
the current-voltage curve over the range 310-345 volts, that is, cor- 
responding to the rounded portion at the foot of Curve D. Curve 1 was 
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plotted two days after drying had begun, that is on the same day as 
Curve B. The currents over that range of voltage were then small so 
that no auxiliary capacity was used. Curves 2 and 3 were plotted after 
16 and 25 days drying respectively. Here and also in the case of Curves 
4 and 5, acondencer of .001 mf was used in parallel with the electrometer. 
For purposes of comparison with the other curves, the ordinates of Curve 1 
should therefore be divided by 63. While the effect of drying as shown 
by Curve 2 has been to preduce a big increase in the current values, this 
increase had not been maintained in Curve 3— in fact there is a consider- 
able reduction. From this and other evidence it was concluded that 
drying was not going on at a satisfactory rate; certainly not as rapidly 
as in previous experiments.* It was thought possible that the difference 
was to be attributed to the glass vessel as metal only had been used in 
the previous work in which a high degree of drying was shown by the 
rapid appearance of ions of high mobility. The problem is perhaps not 
so much to dry the gas as to dry the walls of the vessel. A metal vessel 
was therefore substituted for the glass bell-jar, fresh P2O; was introduced 
and the apparatus sealed up once more. Curves 4 and 5 were obtained 
with the new apparatus 31 and 39 days later. These curves, which are 
very nearly identical, show a considerable increase in current values over 
those obtained with the first apparatus, but even in this case it is believed 
that a really high degree of drying was not reached. It is probable that 
there were minute leaks in the apparatus and that moist air from outside 
found its way in. 

The curves which have been reproduced in Fig. 3 are selected curves 
only in the sense that curves have been picked which cover the whole 
range. Every electrometer observation has been plotted. These curves 
as well as the much greater number which have not been given, agree in 
showing that the current-voltage curve is made up of a number of curves, 
convex upwards, showing intersections at 317, 327 and 337.5 volts. In 
the absence of any explanation of how a periodic character could be im- 
pressed upon the curve by any property of the field, of the grid or of the 
recording instruments, there is no escape from the conclusion that each 
intersection represents a distinct group of ions. 

In Fig. 4 are given six curves covering the range 340-385 volts. As 
before, these curves represent different stages of drying, the lowest curve 
having been plotted 1 day after drying began. All curves are plotted to 
the same scale. Only one of these curves (the short curve on the left-hand 
side) was obtained with the apparatus in the metal enclosure. All curves 
agree in showing discontinuities at 348, 358.5, 368.5, and 378.5 volts. It is 
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hardly necessary to point out that if only a few widely spaced observations 
are taken, the true nature of the curve will not be disclosed. 
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CURRENT 


VOLTS 
540 550 J60 J70 J8O J90 
Fig. 4. Current-voltage curves for 340-390 volts, at different stages of drying. 











From the accurately established critical voltages the following mobility 
values are deduced: 


Critical Voltage Mobility 
317 2.24 
327 2. 
337.5 .10 
348 .04 
358.5 .98 
368.5 .93 
378.5 .87 
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It is not proposed in the present paper to deal in detail with the ex- 
amination of the curves for higher or lower mobilities than those just 
given. The values for the mobilities of the faster ions are not yet all 
definitely known, and what is perhaps more important, the number of 
groups intervening between the fastest given above (2.24) and the ion 
of mobility about 12 (which the author regards as the mono-molecular 
ion) is still slightly uncertain. Judging from the best and most recent 
experiments, the first three ions have mobilities about 11.5, 7.5 and 6 
and there appear to be ten or eleven other ions intervening between these 
and the fastest ion given in the table. As regards the slower ions there 
is ample evidence. Curve B has already shown that there is an inflection 
for about every 10 volts between 250 and 470 volts. Thus the list above 
can be continued giving critical voltages at about 390, 400 etc. to 470 
volts. This is amply confirmed by observations taken over this range 
with both positive and negative ions. Accurate observations show that 
the critical voltages should be approximately 388, 398, 408 etc. and that 
this succession continues up to 528 volts and beyond it. Of course, if the 
air is very dry, Curve C shows that the slower groups do not appear. 

The author interprets these results in terms of a theory of the con- 
stitution of the ion put forward by him some years ago.* According to 
this theory the ion is a cluster of water molecules containing any number 
of molecules up to 30 or 40. Any size of cluster may be present depending 
on the conditions; for example in moist air four types are generally found. 
Drying tends, especially in the case of the negative ion, to eliminate the 
larger clusters and to bring the smaller ones into prominence. In the 
original paper the author suggested that the four kinds of ion ordinarily 
prominent in moist gases contained 18, 24, 30 and 36 molecules of water 
respectively. The mobility of these four hypothetical ions was calculated 
using a formula due to Sir J. J. Thomson. The value found for the 
mobility of the ion composed of 30 water-molecules was 1.53. This ion 
is usually very prominent; Loeb’s Curve E shows it having mobility 1.47 
and Curve A shows it having saturation voltage 474 and mobility 1.50. 
The close analysis of the curves shows that the correct saturation voltage 
should be 468 (mob.=1.52). If our view as to the constitution of the 
ions is correct, the next slower ion—critical voltage 478—has 31 molecules 
of water, the next 32 and so on. If we continue in this fashion we find that 
the ion of 36 water-molecules should be disclosed at critical voltage 528. 
This is a notably prominent ion in all curves taken with moist air. Its 
mobility is 1.34, which corresponds with the calculated value 1.31. Its 


§ J. J. Nolan, Proc. Roy. Soc. 94, 112 (1917) 
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prominence is shown by Curve A, which gives its mobility as 1.36. On 
the other hand, going from the ion of 30 water-molecules to the smaller 
ions, we find that the ion indicated by the critical voltage 408 should have 
24 water-molecules. Critical voltage 408 gives mobility 1.73; the mobility 
calculated by the author for an ion of 24 water-molecules is 1.72. Con- 
tinuing in the same way, the ion indicated by the critical voltage 348 
should have 18 water-molecules. Mobility determined (see table) = 2.04; 
mobility calculated =1.97. The fastest ion in the table should have 15 
water-molecules; the value calculated is 2.17, the value found is 2.24. 
There should then be 14 possible ions left; it has been already stated that 
that number, or possibly one less than it, have been found. 

The numbers (18, 24, 30 and 36) originally suggested for the more 
prominent ions were chosen because they fitted the results of experiment 
better than any others. ‘“‘With fuller information the ions might be better 
represented in terms of some other unit rather than six.’”’® There would 
now seem to be considerable justification for believing that the numbers 
originally suggested cannot be far wrong. Apart from the evidence just 
given, the work on gaseous hydrates® which shows that bodies of the type 
M(H20O), are readily formed, even when M is a molecule of such a gas 
as argon, is highly significant in this connection. 


DISCUSSION 


The author has in previous papers shown by various methods that the 
ions in moist air fall into four distinct groups. Direct confirmation of this 
has been found in the work of C. T. R. Wilson. It is now possible by re- 
fined observation to follow the building up of the ions constituting these 
groups and to detect the less prominent intermediate stages. The num- 
ber of intermediate stages found is in accordance with the number of 
molecules previously assigned to each ion. The smaller ions, when brought 
into prominence by extreme drying, afford further confirmation by their 
absolute mobility values and by their spacing. For many reasons it is 
believed that the unit out of which the ions are built is the water-molecule. 

Apart from the evidence brought forward in this paper, the theory 
which regards the ordinary ion as a charged gas-molecule is quite un- 
tenable. When the author originally found! an ion of mobility 12.5 in 
the air ionized by the decay products of radium, he was unaware of the 
results of Franck,'® Altberg".and Moore.'!* Franck found a negative ion 


®See Bouzat, C. R. 176, 253 (1923), and de Focrand, C. R. 176, 355, (1923) 
1° Franck, Ann. der Physik 21, 972, (1906) 

" Altberg, Ann. der Physik 37, 849 (1912) 

12 Moore, Phys. Rev. 34, 81, (1912) ~ 
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with a mobility of 12.26 in dry air; Altberg found a negative ion of 
mobility 10.1 and Moore found a negative ion of mobility 9.72. The 
agreement is good considering the great differences in the methods em- 
ployed by these workers. That this is the true monomolecular ion can 
hardly be doubted. The values found for the fastest positive ions by these 
workers are respectively 3.22, 7.5 and 6.43. Again, the present writer 
can find no reason for refusing full credit to the results of Haines" on . 
hydrogen; his own observations confirm those of Haines and he has re- 
cently in dry hydrogen found Haines’s ion of mobility 40.6. Finally, the 
fact that the recoil atom has about the same mobility as the ordinary 
positive ion is fatal to the ‘“‘small-ion’”’ view, unless a form of theory is 
invoked which would make the mobility independent of the mass. Even 
this very questionable alternative is ruled out by the observations of 
C. T. R. Wilson® who has obtained photographs showing definitely that 
the atom of ThA was displaced twice as far by an electric field as the 
ordinary ions in the same time. The fact that a body of mass 216 can move 
twice as fast as the ion under an electric field, shows not only that the ion 
is a cluster, but that it is a cluster containing a considerable number of 
molecules. 

The author is indebted to Mr. P. B. Carolan from whose observations 
the curves in this paper were plotted. It is hoped to publish the results 
of Mr. Carolan’s observations more fully on a subsequent occasion, es- 


pecially those dealing with ions of high mobility which have been only 
lightly touched upon in this paper. 
UNIVERSITY COLLEGE, DUBLIN. 
December 8, 1923. 


* Haines, Phil. Mag. 30, 503, (1915) and 31, 339, (1916) 





PRODUCING A SQUARE WAVE OF RADIO FREQUENCY 


A METHOD OF PRODUCING A SQUARE WAVE 
OF RADIO FREQUENCY 


By J. L. Bowman 


ABSTRACT 


A radio-frequency square wave, suitable for the measurement of the mobility 
of ions, has been obtained by using a circuit containing three vacuum tubes. 
The first tube operates as an oscillator. The potential of the second tube, which 
operates as a rectifier, is controlled by the grid circuit coil of the oscillator. 
The third tube is connected so that the plate current will be cut off during 
the half period when the rectified current is flowing in the plate circuit of tube 2. 
The rise of the voltage of the square wave is approximately exponential, due 
to the charge on the grid and the plate of tubes 2 and 3, which must be dis- 
charged through a resistance which couples the plate circuit of tube 2 to the 
grid of tube 3. A cathode ray oscillograph tube indicates that the voltage 
rises to approximately its maximum value in 1/10 of the half period. A square 
wave of 120 volts may be obtained with Western Electric ‘“E’’ tubes. Three 
methods of amplification are given by which a square wave with a very much 
larger voltage may be obtained. 


INTRODUCTION 


HE apparatus for producing a square wave of radio-frequency has 
been designed because of the application it has in the measurement 
of the mobility of ions. The investigations have been confined to vacuum 
tube methods. The circuits found to be most satisfactory are shown in 
Fig. 1. Tube 1 operates as an oscillator, so arranged that the mutual 















































Fig. 1. Diagram of circuit. 


inductance of the coils A, B, C, D, may be varied. Tube 2 may be 
regarded as a rectifier and the circuits immediately associated with tube 3 
produce the square wave. 
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MrETHOD OF OPERATION 


The method of operation is as follows. The oscillator imposes an 
alternating e.m.f. on the grid of tube 2. A certain constant negative 
potential is placed on the grid by the battery /,2 such that during almost 
one-half of the period the alternating e.m.f. causes the plate current to 
be cut off. Variations of the plate current impose an alternating e.m.f. 
on the grid of tube 3. During the time there is no current through R,;, 
the grid of tube 3 is zero with respect to its filament, and a constant cur- 
rent will flow through R,;, but during the other half period there will be 
a large negative potential imposed on the grid of tube 3, thus cutting off 
the plate current. Tubes 2 and 3 must work below the characteristic 
curve during a part of the period. The square wave thus consists of a 
sudden rise to a constant voltage and a rapid decrease to zero voltage. 

The wave will be very nearly a perfect square if the e.m.f. imposed 
on the grid of tube 2 is large. This end may be secured without injury 
to the tube by arranging the circuits as shown in Fig. 1. The grid of tube 
2 connects to the grid of tube 1, the potential of which is controlled by 
the alternating e.m.f. across the grid circuit coil of the oscillator. During 
the negative half period the grid of tube 2 will reach a very low negative 
potential, but during the positive half period the rise will be very rapid 
to about 15 volts, when a current from grid 1 and grid 2 opposes the in- 
crease. Therefore the positive half cycle does not increase the grid 
abnormally. This half cycle, having been amplified and inverted, is im- 
posed on the grid of tube 3. The result is that a square wave of large 
voltage may be obtained without using an amplifier. 


THE FORM OF THE SQUARE WAVE 


The form of the square wave was observed by applying a part of the 
voltage of the square wave to one nair of plates (O, P) of a Western 
Electric cathode ray oscillograph tube (Fig. 1) and by connecting the 
other pair of plates (7, N) to the ends of a resistance in series with a 
variable condenser and an inductance, which is inductively coupled to 
the oscillation circuit of tube 1. By proper adjustment of the resistance 
R’, the capacity C’ and the mutual induction of coils E and F, the e.m.f. 
imposed on one pair of plates of the oscillograph tube will be about 90° 
out of phase with the e.m.f. of the other plates. The sine wave imposed 
on plates \/ and N produces a straight line on the oscillograph and the 
square wave produces a straight line at right angles to it. When the 


difference of phase is about 90° the voltage changes of the square wave 


take place when the e.m.f. of the sine wave is very near a maximum or a 
minimum. Thus the electron beam may be deflected to the right by the 





PRODUCING .L SQUARE WAVE OF RADIO FREQUENCY 3 


sine wave, then vertically upward by the square wave then to the left 
by the sine wave and then vertically downward by the square wave, the 
entire curve being traced once each period. The photographs of the 
screen of the oscillograph show the form of the curves thus obtained. 
Fig. 2 represents 60 volts and Fig. 3 represents 120 volts. 

The vertical lines in the photograph are quite indistinct, indicating 
that the rise and fall of the voltage of the square wave is very rapid. 
The deflection along the horizontal lines represents a simple harmonic 
movement of the electron beam, so that the voltage of the square wave 
changes near the time when the rate of the horizontal deflection is slowest. 
Therefore the curve obtained in this manner does not-represent the exact 
shape of the square wave which would be obtained by plotting the voltage 
against the time /. 























Figs. 2-5. Oscillograph records. 


If RX’ and C’ are removed from the circuit and the plates M, N are 
connected directly to the ends of the coil F, the sine wave voltage and the 
square wave voltage will be in phase, and the vertical deflections will 
take place while the horizontal deflections are most rapid. The curve 
obtained in this manner is shown in Fig. 4. By measuring the horizontal 
displacement which takes place while the beam is being displaced verti- 
cally, the time of rise of the voltage of the square wave was found to be 
about 1,10 of the half period and the time of fall is about 2/3 of the time 


of rise. From these observations approximate calculations may be made 


of the form of the square wave which is represented in Fig. 5. 
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The form of the square wave depends (1) on the amplitude of the 
wave imposed on the grid of tube 2; (2) on the capacities of the various 
circuits and (3) on the voltage amplification factor of tube 2. Using 
Western Electric ““E”’ tubes, with circuit constants to be given later, the 
time required for the original wave to vary over a range sufficient to 
produce the square wave was about 1 20 of the period of oscillation. 

The rate of rise of the voltage of the square wave decreases as the 
capacity of the circuits and tubes is increased. Due to these capacities 
the voltage changes will be approximately exponential, as the charge on 
the plate and the grid is being discharged through the external resistance 
of the plate circuit. The voltage will increase rapidly at first, but as the 
grid potential approaches a limiting value the rate of increase gradually 
becomes less. A slight bend showing this effect can be observed in the 
photographs of the curve (Figs. 2 and 3). As the voltage of the square 
wave decreases, the grid of tube 3 reaches its limiting potential after the 
current in the plate circuit has been cut off. This eliminates a part of 
the capacity effect, so that the time of fall will be less than the time of 
rise of the voltage. 

The capacity in series with R,; consists of the capacity of the plate 
and the filament and of the plate and the grid of tube 2 and the capacity 
of the grid and the filament and of the grid and the plate of tube 3. 
As the grid varies in potential the plate also varies but with the opposite 
potential, so that the total discharge through R,; is of such magnitude 
that the square wave will be appreciably distorted if R,; should be more 
than 2000 ohms when F;,; is 200 volts, but for higher plate voltages it has 
been found empirically that R,; may be approximately equal to 
Ros+ Ro) woin which Ry is the internal impedance of the tube and yo is the 
amplification constant. The amplification factor « of tube 2 will be small, 
as can easily be observed by plotting the voltage across R,; against the 
grid voltage, or it may be obtained approximately from the relation 
H=woR (R+Ro).' The rate of rise of the voltage of the square wave, due 
to an increase of the amplification factor, is opposed by the capacity 
effects of the tube, thus requiring proper adjustment of R,; to obtain 
the best wave form. 


CONDITION FOR Brest OPERATION 


If we desire to produce a square wave for any assigned purpose it is 
necessary to determine (1) the voltage which is to be obtained from the 
square wave; (2) the current that can be taken off without too great 


*Van Der Bijl, Proc. I nst. Radio Eng. 7, 97 (1919) 
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change of voltage; and (3) the resistance which can be used in the circuit 
to which the apparatus is coupled so that the added capacities will not 
distort the wave form. In order to keep the resistance small the im- 
pedance of the tubes should be small. The best operation is found by 
making Rs; about equal to the impedance of the tube. The voltage of 
the square wave obtainable will then be about 1/2 of Fy;._ It may be 
increased to 2/3 of Ey; by increasing the resistance R,;. Or it may be 
increased by making E,3 positive, say 10 volts or even more as a grid 
current through R,; causes no appreciable disturbance. The tube may 
even be operated up to the saturation point of the curve, which is quite 
low when the resistance Ry; is large. The actual square wave voltage 
thus obtained with ‘‘E”’ tubes, was 200 volts by using 375 volts on the 
plate. The voltage F,; may be increased above its normal operating 
value, since there is a large resistance in the grid circuit and because 
the plate current is flowing only half-one of the time and when it is flow- 
ing the voltage between the plate and the filament is decreased to 
(F,—RI). 

After choosing proper values for £4; and Ry3, R,s should be approxi- 
mately equal to (Ry3+Ro)/uo but it is better to use a variable resistance 
and to adjust it to its maximum value while observing the shape of the 
wave. Tube 1 operates best when £;; is about 180 volts, and tube 2 
operates successfully with a plate voltage of 140 volts. Plate batteries 
are shown in Fig. 1 for simplicity, but they may all be replaced by one 
d.c. generator of 400 or 500 volts. This can be done by placing a large 
resistance in the generator circuit and joining the negative end of the 
resistance to the filament of tube 1 and tapping off the desired voltages 
to the plates of tubes 1 and 2. The point at which the circuit from the 
plate of tube 2 connects to the resistance in the generator circuit is directly 
connected to the filament of tube 3, so that the desired plate voltage for 
tube 3 may be obtained by tapping off at a point of higher potential. 
If an oscillograph tube is used to observe the wave form, the generator 
may also be used to place a suitable potential on the anode similar to 
the connections shown in Fig. 1. The adjustment so that each half period 
of the wave is equal, is made by varying the voltage on the battery E,2, 
which should be about 20 volts negative. The oscillator should produce 
a voltage across the condenser of about 450 volts. The frequency used 
was about 100,000 but the wave form is almost as good up to 600,000, 
which is as high as trials were made. 

The resistance should be as non-inductive as possible, though ordinary 
resistance boxes are quite successful. Very good results were obtained 
by using resistances made of quartz fiber sputtered with platinum and 
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mounted in oil.?. It is essential that capacities in ali circuits be made as 
small as possible. If the connections are made as shown in Fig. 1, the 
capacity effects between the batteries will be largely eliminated. If the 
batteries are replaced by a generator, a condenser should be placed in 
parallel with each part of the resistance of the generator circuit which is 
placed in a plate circuit. 


METHODs OF AMPLIFICATION 


If a square wave having an e.m.f. of more than 150 volts is desired, 
amplification may be obtained by either of the three following methods. 
The voltage of the rectified wave of tube 2 may be increased by placing 
another tube in parallel with it, making it possible to increase the voltage 
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Figs. 6-7. Circuits used for amplification. 


on the plate of tube 3. Another method (Fig. 6) is to place the extra tube 
so that the grid is operated from a part of the voltage across R,3, and to 
connect the plate in series with R,;. In this way the negative range of 
potential of the grid of tube 3 is greatly increased, since the current in 
the plate circuit of tube 4 is added to the current of the plate circuit of 
tube 2. In fact the amount of current fed back into the grid circuit of 
tube 3 can easily be adjusted so that a very small range of the original 
wave imposed on the grid of tube 3 is necessary to produce a square wave 
of good form. This method will be very useful if for any reason the 
amplitude of the original wave is small. The voltage of Ey; should be 


? Frayne, Phys. Rev. 21, 348 (1923) 
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about 1/2 of #3. Another method of amplification is to add another tube 
(Fig. 7) to be operated similarly to tube 3 by imposing a part of the voltage 
across Ry; on the grid of tube 4. The total voltage obtained by this 
method, when R,; and Ry; are placed in series in the circuit in which the 
square wave is to be used, will be the sum of the voltages across Ry3 and 
Ry. If this method is combined with the first method of amplification 
given, a very large voltage may be obtained. 

In conclusion I wish to express my thanks to Prof. W. F. G. Swann 
for suggesting the problem and for the valuable assistance he has given. 
I also wish to thank Dr. F. M. Kannenstine for the suggestions he has 
offered. The investigation has been made possible by the courtesy of 
Dr. E. H. Colpitts of the Western Electric Co., who loaned the tubes 
which were used in the experiments. 


RYERSON PHysICAL LABORATORY, 
CHICAGO UNIVERSITY. 
March 1, 1924 
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THE THEORY OF THERMIONICS 
By H. A. Witson 


ABSTRACT 


Thermodynamical theory of electron emission from metals.—By means 
of a reversible cycle similar to that first used by O. W. Richardson, an accurate 
expression for the thermionic emission is obtained and it is shown that the 
constant in the equation for the vapor pressure of the electrons is equal to the 
chemical constant for electron gas. It is shown that if all the internal heat 
of evaporation of the electrons is due to electrical forces then Richardson’s 
original equation J =A 6'/2e~/® must be correct and the specific heat of the 
negative el-ctricity in the metal must be the same as for a monatomic gas at 
constant volume, but that if this specific heat is small, as is known to be the 
case, then all the internal heat of evaporation cannot be due to electrical forces 
and Richardson's second equation J = A 6%«~°’* is a better first approximation to 
the exact relation. It is shown that the equation L = R8@,d(log p)/d@ where L is 
the heat of evaporation of a liquid at pressure p and temperature @ is not exactly 
true for the evaporation of electrons. It is pointed out that the thermo- 
dynamical theory of evaporation was applied to electron emission by the writer 
in 1903 and that O. W. Richardson in 1912, M. V. Laue in 1918, R. C. Tolman 
in 1921 and S. Dushman in 1923 have made further applications of it to the 
same question. 

Quantum theory of the chemical constant.—A simple way of calculating 
the size of the region element in the state space for a monatomic gas is given 
which makes it equal to Nh where N is the number of molecules in the gas 
and fA Planck’s constant. This gives the same expression for the chemical 
constant as that found by Sackur. 


HE term thermionics was introduced by O. W. Richardson to denote 

the branch of physics dealing with the emission of electrons by 
metals at high temperatures. Richardson’s kinetic theory of the emis- 
sion appeared in 1901. 

The emission of electrons by hot bodies is a kind of evaporation and 
the thermodynamical theory of evaporation was applied by the writer! 
in 1903 to calculate the variation of the emission with the temperature. 
The electrons were supposed to behave like a perfect gas which could be 
in equilibrium with the metal at a pressure p and the latent heat of 
evaporation L of one mol was taken to be given by the well known 
equation* 


d (log p) 


L=Ro-——— 


1H. A. Wilson, Phil. Trans. A 202, 243 (1903) 


e . . . . . 
It is shown in the present paper that this equation is not exactly true for electron 
gas evaporating from a metal. 
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where R is the gas constant for one mol and @ the absolute temperature. 
Integration of this equation gave p in terms of @ and the emission per cm? 
per sec. was obtained from p by supposing that the kinetic theory of a 
perfect gas could be applied to the electron gas and that the number of 
electrons emitted was equal to the number striking unit area in the gas 
at pressure p. The internal heat of evaporation P was supposed to be 
given by L=P+Ré@ and P was supposed to vary with the temperature so 
that P= Po+a0, where a is a constant. The pressure p was then given by 


at (S41)! 6+ constant 
og p= Re R og constan 


and the thermionic current J by 
[=zA gita/R .— Po/ RO 


where A is a constant. 

Richardson? considered a reversible cycle in which electrons were sup- 
posed to evaporate from a metal at a high temperature and to be con- 
densed in the same metal at a lower temperature. In this way he got 


l a ane + [+ tant 
og p= Re 2 og Ro constan 


where S is the heat capacity of the electricity in the metal per Faraday. 
This equation, if we take S to be constant, shows that the constant @ in 
the writer’s equation above must be equal to $R+S. On Richardson’s 
original theory we should expect to have S=—%R so that a=0 and 


P 
log p= ating 6+ constant 


which gives Richardson’s well known original equation 


[=A Gi? <~ Po RO 
Richardson pointed out that the values of S found for different metals 
are all small compared with R so that it is approximately correct to put 
S=0 and using also L = R6°d(log p)/d@ he showed that P= P)+$R8@ and 

<> 5 ter et-cnnent 

og p= ——+-lo constan 

hl ek Yee 
which gives 
T= A Gre PFO, 


Richardson definitely adopted this equation as theoretically superior 
to his original equation, though both represent the experimental results 


2 Richardson, Phil. Mag. 23, 602 (1912); ‘Electron Theory of Matter,” p. 448 
Cambridge 1914; also ‘‘Emission of Electrons from Hot Bodies,” p. 30, Longmans 1921- 
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equally well, and he proved that the constant A must have nearly the 
same value for all metals. Later® he found that this equation can be 
obtained by an application of the quantum theory to the problem and 
obtained a theoretical value for the constant A. 

In a recent paper on electron emission from metals S. Dushman‘ has 
calculated the emission thermodynamically by assuming it analogous to 
evaporation and has determined the absolute value of the constant in 
the equation for log p by putting it equal to the chemical constant for 
electron gas as calculated by Sackur, Tetrode and Stern by means of 
Boltzmann and Planck’s theory of entropy. Dushman also assumes that 
the heat capacity of the electrons in the metal is negligible and so gets 
Richardson’s second equation J = A 6? e~?*/®? with A =2xmek?/h® where k 
is the gas constant for one molecule, m the mass and e the charge of an 
electron and h Planck’s constant. 

The assumption that the thermodynamical theory of evaporation can 
be applied without modification to thermionic emission requires justifica- 
tion and Richardson’s reversible cycle supplies this up to a certain point 
but it is not clear that the constant in the equation for the vapor pressure 
of the electron gas is equal to the chemical constant for electrons. Tol- 
man® has shown that the entropy of the electron gas calculated from the 
quantum theory of entropy agrees nearly with that calculated from the 
heat of evaporation, which shows that Dushman’s assumption is at any 
rate approximately true. 

The quantum theories of the chemical constant are not satisfactory. 
The main objects of the present paper are to attempt to show that the 
constant in the equation for log p is exactly equal to the chemical con- 
stant for electron gas, to give a simple method of calculating the chemical 
constant on the quantum theory and to obtain an accurate expression 
for the thermionic emission on the thermodynamical theory of it. 

The analogy between the evaporation of a liquid or solid and the emis- 
sion of electrons by a metal is not perfect. If any appreciable number of 
electrons is supposed emitted the negative charge in the gas and the 
positive charge on the metal become too large for their effects to be 
neglected. This difficulty can be overcome by using a reversible cycle 
in which electrons are supposed to be emitted at one temperature and 
simultaneously absorbed at another temperature so that the number of 
electrons outside the metal can be supposed constant and very small. 
A cycle of this kind was used by Richardson but he did not obtain the 

* Richardson, Phil. Mag. 28, 633 (1914) 


* Dushman, Phys. Rev. 21, 623, (1923) 
® Tolman, J. Bur. Chem. Soc. 43, 1592 (1921) 
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complete solution of the problem by means of it alone. The application 
of the first and second laws of thermodynamics to such a cycle is sufficient 
for the solution of the problem but the equations obtained differ from 
those for the evaporation of a liquid by terms containing the electrical 
potentials. By supposing these potentials equal to zero the equations 
become precisely similar to the equations for evaporation of a liquid so 
that the constant in the vapor pressure equation can be shown to be 
equal to the chemical constant for electron gas. 

It will be convenient to compare two similar reversible cycles one in- 
volving the evaporation and condensation of a liquid and the other 
the evaporation and condensation of electrons from a metal. 

Consider two vessels containing the liquid one maintained at constant 
temperature 6, and the other at 6, and let them be connected by a fine 
tube through which the liquid can flow. Let the liquid at 6, be in equi- 
librium with its vapor at pressure p, and let that at 62 be in equilibrium 
with its vapor at 2. Now suppose one mol of liquid evaporated from 
that at 6, and one mol of vapor condensed in the liquid at 62 and let one 
mol of the liquid flow through the tube from 6, to 6. The heat required 
to maintain the temperatures constant during this operation will be 
L,—L2.+fcd0( 62 to 0;) where L; and L» are the latent heats per mol at 
6, and 6, and cis the molal heat capacity of the liquid.. The system does 
work R(6,— 62) if we assume the vapors can be regarded as perfect gases. 

Now cool one mol of the vapor at 6, and ~; down to 6 at constant pres- 
sure ~, and then expand it at constant temperature 6, to pressure p» 
and add it to the vapor in equilibrium with the liquid at #. The heat 
required for this operation is Cp(02—6:)+R62 log (p;/p2). The work done 
by the system is R(6:— 6;)+ R62 log (p:/p2). The final state of the system 
is the same as at the start so that we have 


61 
Li-lat f cdé — R(0, — 62) +C,(62—6;) —R(82—6;) = 0. 
Os 
Hence 


a 
La=LatC(-6)— f cae 
6: 


Also the total change of entropy is zero so that 


7 I, ln 6: cdO 
conatennaiie ——C,Jog— 
6; Oe 6: 6 6 
If 6.=0 equation (1) gives, on putting C,=§R which is true for mon- 
atomic gases, 


6 
“+R log Pig (2) 
r p 


5 A 
Lelet; Ré - [can (3) 
0 
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and from (2) and (3) we get 


Lo a cd 1 
oR log 4+R log pi+ a ee cd6= 


1 6; 0 


Lo 9s cdé td 
Le Seg et 62+ R log po+ eee server 


2 


a a. f ca0- fs (4) 
og p= Re 7 °8 C 


where # is the constant usually called the chemical constant for the liquid. 

Now consider a similar reversible cycle for electrons emitted by a 
metal. Let there be two pieces of the metal maintained at temperatures 
6, and 6. and connected by a fine wire of the same metal. Let them be 
in equilibrium with electron gas at pressures p; and p2. We suppose these 
pressures so small that the mutual repulsion of the electrons can be neg- 
lected and the electrons regarded as a perfect monatomic gas. Suppose 
a reversible cycle of operations carried out with the electron gas, pre- 
cisely similar to that described above. Let the electrical potential inside 
the metal be zero and let it be ¢; just outside the metal at 6, and ¢» just 
outside that at @. The latent heat of the electrons will be taken to be the 
heat absorbed at the surface of the metal when one mol of electrons 
escapes, the potential of the electrons changing from zero in the metal 
to ¢ outside and the potential of the metal being kept at zero by allow- 
ing electrons to flow into it through the wire. We suppose all the opera- 
tions of the cycle carried out simultaneously so that the number of 
electrons outside the metal remains constant during the cycle. Then 
when one mol of electron gas is transferred from the gas at 6, to that 
at 6. work Ne (¢2:—¢,) has to be done against electrical forces, where e 
is the charge on one electron and N the number of electrons in one mol. 
Thus we have 


5 6 
L=Lot Nelo) +> R0+ f Sd6 (5) 


3 8 
or Pervtine~tt- Rot fas (5a) 
: 0 


where S is the heat capacity per Faraday (Ne) of electricity in the metal 
and P is the internal heat of mens Also as before 


Ilo Ne(o— bo) s, 
] = —— — —____ g d—-— Sd0+ ae 
08 P Re Re ak ’ i , (6) 


where 7’ is a constant. 
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Equation (6) differs from the corresponding equation (4) by the term 
— Ne (¢—¢0)/R@ but the constant 7’ can nevertheless easily be seen to 
be equal to the chemical constant for electron gas. 

The latent heat of evaporation of the electrons is partly due to the po- 
tential difference between a point just outside and one inside the metal. 
We may suppose L = w+ Ne@ where w is the part of the heat of evapora- 
tion due to other than electrical forces. Equation (5) then gives 


5 6 
w=ut= Rot f Sao (7) 
0 
and (6) gives 


| mn wt hog + 1 [saote. 8) 
8 P= — Ro Re RoJo 7 


The change of potential ¢@ at the surface of the metal may be supposed 
to take place in a double layer at the surface so that if p’ is the pressure 
just inside this layer we have 

log p’=log p+ Ned/RO 
Hence 


0Sd0 


@ 
] = —— a 0 ——— | Sd6+7’. 
06 P et ” +f R6@ R@ ie 


This equation gives the relation between the temperature and the vapor 
pressure which the electrons would have if they were uncharged or if 
@ were zero and wp is the heat of evaporation they would have at the 
absolute zero if @ were zero. 

This equation is therefore precisely similar to equation (4) which gives 
the vapor pressure of the liquid, so that it is clear that the constants 7, 
and 7’ are equal respectively to the chemical constants of the liquid and 
of electron gas. Hence we have 

(Lo+Ne(o— do) +fpSd0) 5 fz 4 


l =— —log 0 — 9 
og p Re +; og 0+ +1 (9) 


where 7 is now the chemical constant for electron gas. 
Equations (9) and (5) give 





L Ren )+N ot 
= (8) =. 
do 8? d6 


Thus it appears that the equation 


d 
L=Re—(I 
70 \°8 p) 


which gives the latent heat of evaporation of a liquid is not exactly true 
for the evaporation of electron gas. The previous thermodynamical 
theories of the emission of electrons given by the writer, O. W. Richardson 
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and others in which this equation was used are therefore approximately 
correct only when the term Neéd¢/d@ is relatively small. 

In Richardson’s® original theory it was supposed that the internal 
heat of evaporation was entirely due to electrical forces so that wo=0 
and w=R8. In this case equation (7) gives 


3 @ 
—- R@= - f sas 
2 0 


so that S=—3R. This is the”value off that we should expect if the 
electrons in the metal had the same energy as gas molecules. Putting 
this value of S in (6) or (9) we get 

log p= — Ned/RO+ log 6+ constant (10) 
which gives Richardson’s original equation for J. But as Richardson 
pointed out, the values found experimentally for S are small compared 
with R so that as a first approximation we may put S=0 and then (7) 
gives w=wo+$R0 so that w cannot be less than $R@ and (9) becomes 


Io+Ne(o— $0) 5 
- —log 6+72. 11 
os +; og 0+1 (11) 


If we also suppose that ¢ is independent of the temperature so that 
$— do=0 we get 





log p= — 


] = a 6+2 (12) 
og p= Ro 2 8 1 


which is the equation used by Dushman. Lp» of course is equal to Po or 
to w+ Ned¢o. The theory of thermo-electricity indicates that Sis probably 
proportional to the absolute temperature @. If we assume S=.S,6 where 
Si is a constant, equation (5) becomes 


5 1 
slat NAO— te) +, RO+ Se (13) 


and (9) gives 
Lot+Ne(o— do) 5 S, 6 
-- — log 6+—+7. 14 
a +; og tet (14) 


If J denotes the thermionic emission current we have 
pe(1—r) 
2xmké 


where k is the gas constant for one molecule and r is the reflecting power 
of the metal for electrons. This equation with (9) or (14) gives J in terms 
of 6. 





log p= 


I (15) 


* Richardson, Camb. Phil. Proc. 11, 286 (1901) 
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The chemical constant 7 was calculated on a quantum theory by Sackur’ 
and found equal to i9+$log M where M is the atomic weight of the 
monatomic gas. Tetrode*, Stern’, Planck'® and Brody" have confirmed 
this result but get a different value for 7%. Lewis'* favors a value of ig 
which lies between those of Sackur and Stern. All these values are claimed 
to be supported by the experimental evidence. Since the theory relates 
to an ideal perfect gas whereas the experimental evidence is got from 
experiments on real gases which are not perfect we should not expect 
exact agreement between the theoretical and experimental values. It 
seems likely that the experimental values should be a little greater than 
the theoretical values for perfect monatomic gases. If this is admitted 
then Sackur’s theoretical value is most probably the correct one for a 
perfect gas because it is a little smaller than the experimental values 
whereas Stern’s value is rather larger. 

The assumptions made in the calculations of 7 just mentioned do not 
inspire much confidence so the following simple method of calculating 7, 
which is based on the established principles of the quantum theory, may 
be useful. The result obtained agrees with that of Sackur. 

Planck in his ‘“Theory of Heat Radiation’ Chapter II, calculates the 
entropy of an ideal monatomic gas by means of his and Boltzmann’s 
probability theory of entropy and shows that 


: Nk*!2(22m)3/? 
ilog{ } 





eG 


where G is the size of the region element of the state space for N mole- 
cules. It follows as Planck points out that G is proportional to N so that 
if g is the size of the region element for one molecule then G=Ng. This 
result obtained by Planck before 1911 agrees with Sackur’s later result 
if we put g=h*. Planck however did not then determine the value of 
g and did not show that it was the same for any monatomic gas. 

To determine g imagine a rectangular box with smooth rigid walls 
containing a single molecule of a monatomic gas. Let the sides of the box 
be of lengths a, b, and c. Take the origin of coordinates (x, y, z) at one 
corner of the box and let the sides a, b, c lie along the axes x, y, 2 respec- 
tively. The molecule will move in the box along straight lines and when 


7 Sackur, Ann. der Physik 40, 67 (1913) 

§ Tetrode, Ann. der Physik 38, 434; 39, 255 (1912) 

* Stern, Zeit. f. Electrochemie 25, 66 (1919) 

10 Planck, Berlin Akad. 1916 p. 653 

" Brody, Zeit. f. Physik 6, 79 (1921); cf. Bull. Nat. Research Council, “Quantum 
Theory,’’ Adams 

® Lewis, Phys. Rev. 3, 92 (1914) 
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it strikes a side of the box its velocity component perpendicular to that 
side will be reversed. The state space for this molecule will be of six 
dimensions mx, my, mz, x, y, and z. The path of the molecule in the 
(mx, x) plane goes around a rectangle ABCD, with AB and CD parallel to 
the axis of x. At A the molecule is at the end of the box where x=0 
and then moves along AB since x is constant. When x =a it collides with 
the other end of the box at B and its velocity « is reversed at P so the 
path goes to C. The path then goes along CD and back to A when the 
collision at x=0 takes place. The motion is thus a cyclic one and the 
frequency is x/2a. According to the W. Wilson-Sommerfeld quantum 
condition which is supposed to apply to any cyclic motion of a single 
electron or atom we have 

S pdq=S midx= nh 
where m is an integer. The area ABCD must therefore be mh so that 

2mxa = mh 

or x =mh/2ma. 
In the same way we get y= %2h/2mb and z=ns3h/2mc. The kinetic energy 
of the molecule is therefore 


e F ¢ 

This result was obtained by Brody" by a similar method. The region 
element on the plane (mx, x) is therefore of area h and it is made up of 
two rectangular strips of length a and total breadth h/a. The six dimen- 
sional region element is therefore a block with sides a, 6, c, h/a, h/b and 
h/c so that its volume is h*. Hence G= Nh’ and 


: {x an 
t= log : 


e ° - h? n,2 nN” ns 
3m (x?-+y?-+2") -—_( 
8m 





eh’ 
This result becomes the same as Stern’s which was adopted by Dushman, 
if the € is omitted, so it is less than Stern’s value by unity. 

Brody’s calculation which gives Stern’s value of ip really shows that 
Stern’s value is too large. He takes the number of ways in which a gas 
can have energy exactly E equal to the number in which its energy can 
be between E and E+1 which clearly enormously overestimates the num- 
ber of ways and so must give too large a value of the entropy and hence 
also of io. 

Using the value just found for 7 we obtain from equations (9) and 


(15) the exact equation for the thermionic emission current J, 


eoaee-cogieae fz 

2xmek?(1—r) ., - Ro ~J ie (16) 

= Ge . 
eh’ 
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If we suppose the metal does not reflect any electrons so that r=0 and 
that ¢—¢0=0 and also that S=0 we get 
2xme k? 


eh’ 
If Stern’s value of 7 is preferred the ¢ in the denominator may be omitted 
from (16) and (17). 

The reflecting power r is probably not large for slow electrons and ¢@ 
probably does not vary much with the temperature so that this investiga- 
tion may be said to confirm (17) as a first approximation for J. 

The very accurate measurements of J which have recently been made 
by Davisson and Germer," Schlichter™ and others seem to justify the 
working out of an accurate expression for it. 

Dushman shows that the results obtained for tungsten and platinum 
when put in the equation 


«Sel (17) 


2xme k? 
specail 


[= rr \ /0 (18) 


give very constant values of bo over considerable ranges of temperature. 
If we suppose that 


8 0Sdé 
Neo—en)+ [sao—0 f so {1-tog¢t—1)} =0. (19) 
0 


Then since b)>=Lo/R (16) becomes identical with (18) so that if (19) is 
true then Dushman’s equation is exact. Davisson’s results for tungsten 
with equation (18) give values by) which do not vary from the average 
value 51860 by more than 60 for temperatures from 1935 to 2306. The 
value of the specific heat of electricity for most metals is about — 1000 
ergs per e.m.u. at ordinary temperatures. This gives for S about —0.2 


calories per Faraday. If we suppose S =.5,@ then S; must be about —0.001. 
2000 2300 6 


This gives [Sao = — 2000 calories and, /Sdo = — 800. 0f is about double 
0 0 


1900 


4 
S Sao so that the terms containing S in equation (19) are about equal to 


0 
2000 calories and the values of these terms increases by about 800 in the 


range from 1900 to 2300. This is about seven times the variation in 
Lo = 2bo got from Davisson’s results so that to explain them it is necessary 
to suppose that the changes in ¢ and perhaps r with the temperature 
compensate for the changes in the integrals of S. Of course if S is not 
propertional to @ these variations would be smaller but it is clear that 
Davisson’s results are sufficiently accurate to make it worth while to use 


8 Germer, Phys. Rev. 20, 300 (1922) 
“4 Schlichter, Ann. der Phys. 47, 573 (1915) 
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the exact expression for J in place of the first approximation. It seems 
probable that ¢ changes with the temperature in such a way as to more 
or less completely compensate for the changes in the integrals of S and 
the term R@(1—log (1—r)). To test the exact equation (16) it will be 
necessary to have accurate measurements of J over a large range of 


temperature and approximate values of S, r and @¢—¢o for the same 
range. 

¢— do is equal to the contact difference of potential between the hot 
metal and the same metal at the absolute zero. Unless some idea of the 
value of this quantity can be obtained independently it will always be 
possible to make any experimental results for J fit the equation (16) by 
supposing ¢—¢ to vary with the temperature in the necessary way. 
If it is found that equation (18) exactly represents the experimental 
results then it will be necessary to conclude that the relation (19) between 
¢, S,r and @ is true and it may be possible to work out a theoretical basis 
for this equation. 


Rice INSTITUTE, 
Houston, TEXAS. 
January 29, 1924. 
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CURRENTS LIMITED BY SPACE CHARGE BETWEEN 
CONCENTRIC SPHERES 


By IrviNG LANGMUIR AND KATHARINE B. BLODGETT 


ABSTRACT 


Limiting current between concentric spheres; calculation of the function 
a=f(r/ro) in the space charge equation i = (4+/2/9)+/(e/m) V3/2/a?.—The co- 
efficients of the first six terms of a series for a were determined, and a? cal- 
culated from this series. The results were checked by an integration method 
which was also used to calculate values in the region where the series failed. 
For an emitter of radius ro inside a collector of radius r, values of a? when 
log (r/ro) >6.4 are given by the equation 

4a?=0.112 log (log r/ro) +4 log (r/ro) +0.152. 
Where the collector is the inside sphere, values of a? for ro/r>9 are given 
by the equation ($a*)?/?=1.11 (r9/r)—1.64. It is shown that when the col- 
lector is the inside sphere the potential distribution near the collector is un- 
altered if the emitter is replaced by a non-emitting sphere with a diameter 
.677 times the original diameter. 
Limiting current between coaxial cylinders and between concentric 
_ spheres.—Equations are derived for the current in terms of the radius of cur- 
vature of the emitter. It is shown that at a surface in space four-fifths of the 
distance from the emitter to the collector the current density is independent 
of the radius of curvature when r/ro or ro/r<2; and in the case of coaxial 
cylinders with the emitter inside this holds true even when r/r)=20. 


HE problem of the calculation of thermionic currents limited by space 

charge has to deal with three simple cases, those of parallel planes, 
coaxial cylinders, and concentric spheres. These three cases have the 
characteristic that when an electron current is flowing, the lines of force 
and the paths of the electrons coincide. Parallel planes' and coaxial 
cylinders? have been considered in previous papers, and the equations 
have been given for calculating the currents for these surfaces. Although 
it is easy to construct apparatus in which the conditions approximate 
closely those of the cylindrical case, it is difficult to build devices in which 
the heated cathode takes the form of a sphere. 

Recently space charge equations have been used in a new line of ex- 
perimental work* to. measure the intensity of ionization in gaseous dis- 
charges in which a positive ion sheath forms around the negatively charged 
electrode and the current is limited by space charge. Surfaces of different 
shapes have been used for electrodes, and since the sphere has proved 

' Langmuir, Phys. Rev. 2, 450 (1913); Phys. Zeits. 15, 348, (1914) 

? Langmuir and Blodgett, Phys. Rev. 22, 347 (1923) 


* Langmuir, Science, 58, 290 (1923); Gen. Elec. Rev. 26, 731 (1923); Journ. Franklin 
Inst. 196, 751 (1923) 
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to be particularly suitable, the space charge equations for spherical elec- 
trodes have assumed practical importance. The development of these 
equations and the calculation of the numerical values of the function a 
which they involve will be given in this paper. The derivation has been 
made with the intention of using the equations to calculate the flow of 
either electrons or positive ions, and for this reason the general case 
has been developed in which the current flows from an emitter to a 
collector without regard to the direction of the voltage between the 
electrodes. Thus for calculations on ionization in gaseous discharges the 
outer edge of the positive ion sheath itself is the emitter of ions and the 
negative electrode within the sheath is the collector. 

For the case of space charge between concentric spheres Poisson’s 


equation becomes 

$ £2 re 
r? dr r P 

where V is the potential at a point distant r from the center, and p is the 

electron space charge. If 7 is the total thermionic current, v the velocity 

of electrons or ions at the distance r, and e and m the charge and mass of 


an electron or ion, v : 
i=4nr'pV. : (2) 


If we neglect the initial velocities of the particles we have 

3mv?= Ve (3) 
where V is measured from the surface of the emitter. Combining Eqs. 
(1), (2), and (3) gives 


df dV\ . f/m 
— rn) =is/ (4) 
dr dr 2eV 


This equation, like the similar equation for coaxial cylinders, probably 
cannot be directly integrated, but a result can be obtained in terms of a 
series. The final solution takes the form 


472 Se vil | 
nt fe ° 
9 m a? 
where a is a function of the ratio of the radii r and ro of the spheres, ro 
being the radius of the emitter. It is seen from this equation that the 


current is independent of the actual sizes of the spheres since the radii 
appear only in a ratio. 


CALCULATION OF a 


An equation in a can be obtained by substituting Eq. (5) in Eq. (4) 
and placing 
y= log(r/r) (6) 
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This gives 
sect (s +3e2—1=0 (7) 


From this equation the values of a can be obtained in terms of a series. 
The coefficients of the first six terms of the series were obtained by a 
method described in a recent paper,? which makes use of Maclaurin’s 
series in the neighborhood of y=0. The series for a thus calculated is 


a=y—0.3y?+0 .075y'—0 .0143182y'+0 .0021609y' —0 .000267917°+ ...(8) 


The first three coefficients are exact and the last three are each rounded 
off in the last figure. 

Case where r/ro>1 (Emitter inside of collector). For values of r/ro up to 
5, a could be calculated accurately to three places of decimals from the 
series in Eq. (8). For the larger values, however, the six terms were en- 
tirely insufficient, and an integration process was used to complete the 
calculations. This process was applied to the equation obtained from 
Eq. (7) by making the substitution 

y= }a° (9) 

which gives 


(7 *) + +37 -1- 0. (10) 


The method employed was described in the recent paper referred to 
above, and consists essentially in substituting approximate values of y 
and dy/dy in Eq. (10), solving the equation for d*y/dy?, and integrating 
the values thus obtained for new values of y and dy/dy which will be 
more nearly correct than the approximate values substituted at the be- 
ginning of the process. These calculations were made with y in the 
region y = 1.0 to y=2.2, and the values of y obtained as the result of in- 
tegration agreed exactly with those obtained from the series. The process 
was then used to calculate y in the range y = 2.2 to y = 26.8 where the series 
was insufficient. The results of these calculations are given in Table I. 

The following empirical equation was found to fit the calculated values 
of y, dy/dy and d*y/dy? with a high degree of accuracy for values of y 
greater than 11.8, 

y=0.112 log y+4y7+0.152. (11) 


The values of y given by this equation have an error of .48 per cent at 
y= 6.4 and .025 per cent at y= 11.8. 

Table II gives the complete range of values of a? up to r/ro= 100,000. 
The values corresponding to r/ro>5 were obtained from interpolations 
for a from Table I, using Newton’s interpolation formula. 
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Case where r/ro<1 (Emitter outside of collector). The series in Eq. (8) 
that was used to calculate a in the preceding case is also applicable in the 
inverted case where the cathode is the outside cylinder, but since y is 
now negative the signs of the alternate terms are changed. In this case, 
also, the integration process was used to calculate a for the larger values 
of y. 


TABLE I 








dy/dy a 


.4512 . 2980 
.4435 3652 
-4354 4281 
.4275 .4873 
.4198 .5432 


.4126 . 5963 
.4060 .6468 
.3970 . 7183 
3891 . 7856 
. 3825 . 8493 


.3767 .9099 
.3720 .9678 
. 3678 .0234 
.3615 . 1288 
. 3569 .2276 


. 3536 .3214 
.3510 .4107 
. 3493 .4963 
. 3476 .5787 
. 3433 .9532 


.3410 . 2823 
3397 .5799 
. 3387 . 8536 
. 3380 . 1086 
.3375 . 3482 
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The equation which was used in this process was obtained by combin- 
ing Eqs. (6) and (10) and placing ro/r =o, which gives 
d*y 1/fdy\? 1 
3 =--(Z) ——=0 (12) 
do? y\do o 
The results of calculations using the integration process in the range 
o=7 to o=22 are given in Table III. 
In the neighborhood of «= 22, the term 1/e? in Eq. (12) is decreasing 
rapidly and has already become only 0.7 per cent of the second term, 
so that a solution of the equation neglecting the term 1/o? gives a good 
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approximation formula for y which may be used to obtain values of a 
in the region beyond «=22. The complete solution of 


d*y 1 fdy\? 
s—t- (7) =0 (13) 
do? y\do 


y8=Aot+B (14) 
TABLE IIT 
a? as function of radius. 
ro=radius of emitter; r=radius at any point P; 


a? applies to case where P is outside emitter, r>ro; 
(—a)? applies to case where P is inside emitter, ro>r. 
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The values of A and B were derived from the data in Table III, with 


the result 
y8=1.110—1.64. (15) 


Values of y obtained from this equation have an error of 0.5 per cent 
at o=9 and 0.025 per cent at c=18. Table II gives the complete range 
of values of a? for this case up to 7o/r=500. .Eq. (15) was used for the 
calculations in the range 7/r=22 to 500. 


TABLE III 


dy/do y=4a? a 

.0454 .675 3.918 
. 1376 721 4.176 
.2261 .812 4.430 
3113 .947 4.679 
. 3935 .123 4.924 
.4730 3.340 5.165 
.5501 .596 5.403 
.8379 .988 6.323 
0990 .930 . 7.201 
3400 371 8.046 
. 5646 39.280 8.863 
.7760 .622 9.656 
.9761 54.378 10.429 
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Low values of a? for either the direct or inverted case where r/ro or 
ro/r<1.4, may be conveniently calculated from Eq. (16) which is derived 
in the following manner. Put e=(7/r)—1 and expand y = —log(ro/r) 
into a series in terms of e. Substitute this value of y in the equation ob- 
tained by squaring Eq. (8) and as a result 


77 
, (16) 


Logarithmic differentiation of Eq. (5), considering 7 as constant, gives 
da _3 dV 
a 2 Le 


Logarithmic differentiation of Eq. (15) combined with Eq. (9) gives 
4 da 1.11 do 


= aii , (18) 
3 a 1.11l0—1.64 


Combining this with Eq. (17) and substituting ¢ =1r9/r, 


=~ -2/e 1.477 - 


(17) 





In the case of two concentric spheres, neither - which is emitting current, 
the voltage at any point of radius r is 


1 
V=A (---) (20) 
r r; 
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where V is measured with respect to the sphere of radius r;. Logarithmic 
differentiation of Eq. (20) gives 


dV dr (- -) ' 21) 
Vr r nj ( 


Comparing Eqs. (19) and (21) we see that they become identical for all 
values of r if we place 

r,=0.677 ro. (22) 
Thus when current emitted by a spherical electrode is being collected by 
an inner concentric sphere of small diameter, the effect of space charge 
is to make the potential distribution near the collector the same as it 
would be if the outer sphere were not emitting current and were given a 
diameter 67.7 per cent of the original diameter. A similar ratio of di- 
ameters for the case of cylinders was given in the previous paper and is 
0.707. 


COMPARISON OF EQUATIONS FOR PLANES, CYLINDERS AND SPHERES 


The equations- given in this paper and two preceding papers for 
_ calculating thermionic currents limited by space charge are as follows: 
Parallel planes (i,=current per unit area) 


2 ae y3/2 
mL en 
x 2 


Coaxial cylinders (1;=current per unit length) 
y3/2 


rt 


p3/2 ° 
7 (25) 
a? 
If these three equations are written in the same form by substituting 
1,/(2mro) =i, and 1/(42ro) =i,, they become 
for parallel planes, 
i,=D/x?; (26) 
for coaxial cylinders, 
ta=D/(rorB*); (27) 
for concentric spheres, 
tg=D/(r0°a*) ; (28) 


. si 
pew Vertis (29) 
Or m 


where 
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If 7 is expressed in amperes per cm?, V in volts, x and r in cm, and e 
and mare the charge and mass respectively of an electron, D has the value 
D=2.336X10-*V3/? 


If e and m are the charge and mass of an ion, 
D=5.455X10-°V*/2/./M 


where M is the molecular weight of the ions (oxygen atom = 16). 

Consider the case of two concentric spheres that are Very close together 
so that their surfaces are almost parallel. If we put jr=(ro+x) and 
b=x/ro=(r/ro—1) we can write a series for a? similar to the series in 
Eq. (16) as follows 

ata pt pelt bi— (30) 
5 300 © 
This value of a? may be substituted in Eq. (28) and it is then seen that 
when b becomes so small that all the terms in the series except the first 
are negligible, the equation becomes identical with Eq. (26) for parallel 
planes. 

In the same way a series for 8* can be derived, which is 

gt=b?—— 4 b4—... . (31) 
5 50 
and this can be substituted in Eq. (27) which will also become identical 
with the equation for parallel planes when bd is very small. 

Where the surfaces of cylinders and spheres are farther apart so that 
the second and third terms of the series in Eqs. (30) and (31) become 
important and the conditions no longer closely approximate those of 
parallel planes, these series enable us to obtain equations for the current 
in terms of the total curvature of the emitting surfaces. This total 
curvature po in the case of the emitting cylinder is pp=1/7o, and for the 
sphere is pp =2/ro. Combining r=(7o+x) with Eq. (31) and substituting 
the value of po for a cylinder, we have 





rorB? = (1 oa. . :) (32) 
5 50 
Similarly from Eq. (30) we have for spheres 
ra? = aC ea and x* por—.. ) (33) 
5 1200 


We now see that the equations for cylinders and spheres in terms of the 
total curvature of the emitting surface have the first two terms identical. 
By expanding the reciprocal of the series in Eq. (32) we obtain an ex- 
pression for 1/(ror8?), in which the second term is +4/5 (xpo) but the 
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third term is only —0.02 x*p? which is 3 per cent of the third term in 
Eq. (32). Carrying through a similar process for Eq. (33) the third term 
becomes 0.126 x?p*. Substituting these values for 1/(ror8*) and 1/(ro%a?) 
in Eqs. (27) and (28) and neglecting the third terms of the expansions, 
we have 


2 "e ye 4 
I ) (34) 


“ee Va ets 
which is the equation for parallel planes with a first order correction term 
for curved surfaces in terms of the curvature of the surface. 

But it was arbitrary when defining the current per unit area and the 
total curvature to choose the area and curvature of the surface of the 
emitter. We may just as well choose the surface of the collector or any 
surface between them. We shall consider the general case where the 
current is measured per unit area of a surface of radius 7;=(r9+x) for 
which the total curvature is p; = po/(1+ pokx) fora cytnter OF pi = po/(1+ 
3 pokx) for a sphere. 

In this case from Eqs. (27) and (28) the current per unit area of the 
new surface is iz=D/[(ro+kx)r8?] for cylinders and i,=D/[(ro+kx)?a?] 
for spheres. Eqs. (32) and (33) will then become 


(ro+ kx) rp? = x? [1+(#-z)amt( a= be x* pr? ---| (35) 
and 


617 
(rot+kx)?a* = x* e[ie(et)ane(2h e+) stot -—-] . (36) 


It now becomes apparent that for k=4/5 the second term in both 
series will vanish for all values of xp;, and the coefficients of x?p,? will be 
small, having the value .02 for cylinders and .0342 for spheres. This 
means that for all values of xp; that make the third terms negligible, the 
current per unit area of a surface four-fifths of the distance from the emitter 
to the collector is the same as for parallel plane electrodes spaced the same 
distance apart as the two curved electrodes. 

We have thus determined a surface in space at which the current per 
unit area is the same whether the electrodes are plane or curved, provided 
that .03x*p,? is small compared to unity. This suggests that it may be 
possible to locate such a surface in every instance by properly choosing 
the value of &. At this surface, the equations for parallel planes, coaxial 
cylinders, and concentric spheres must become identical, and therefore, 
substituting (ro+x) for ro, we have from Eqs. (26) and (27) 


= (rot kx)rp? (37) 
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from which 
(38) 


Similarly in the case of spheres 
x? = (ro+ kx)?a? (39) 
from which 


(40) 


Since x=(r—7), Eqs. (38) and (40) may be written with the radii of 
the curved surfaces appearing in each case only as a ratio. Values of k 
calculated from these equations are given in Table IV. 


TABLE IV 
Values of k as a function of the radit 








! or 7% = Cylinder Cylinder Sphere Sphere 
ro r/ro>1 r/to<l r/ro>1 r/ro<l 


r 

0 0.8000 0.8000 0.8000 0.8000 
2 .7970 . 8039 . 7877 .8122 
5 . 7932 . 8093 .7718 . 8268 
0 
0 





1 

5. 

B. 

2. . 7904 .8171 .7513 . 8450 

5. . 7935 8455 . 6862 . 8959 
10 . 8090 .8677 .6393 .9260 
20 . 8340 . 8883 . 5959 .9490 

100 .9081 .9258 . 5132 .9805 
oo 1.0000 1.0000 .0000 1.0000 








Consider the case of two coaxial cylinders where the emitter is the 
inner one and has a radius 1 while the collector has a radius 20. From 
Table IV k=.834, and a surface 83.4 per cent of the distance from the 
emitter has a radius r; = 16.846. The current per unit area of this surface 
by Eq. (27) is i7z= D/(16.84678?). Substituting the value 6? = 1.0715 when 
r/ro=20 we have 1,=D/361 for the correct value of the current, which 
is the same value as 1,=D/(19)* for the current between two parallel 
planes. 

If we had given k the approximate value 0.8 we should have had 
1, = D/347.2 for the current per unit area of a surface with radius r; = 16.2. 
This differs from the value of current between parallel planes by 4 per 
cent, and since this error is (.8x —.834x)/(ro+.834x) it is clear that it is 
less than the error in'k. Thus in this case for a range of ratio of radii 
from 1/1 to 1/20 the current per unit area of a surface four-fifths of the 
distance from the emitter is the same as the current between parallel 
planes within an error of 4 per cent, and over half the range the error 
is less than 1 per cent. 
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On the other hand in the case of coaxial cylinders with the emitter 
outside, Table IV shows that k departs much more rapidly from the value 
0.8. The error introduced into calculations of current by choosing k =0.8 
instead of the exact value k, from the Table is (.8x —k,x)/(ro+k,x) where 
x=(r—ro) is now negative. This error in current is proportional to the 
error in k as ka is to (%9/x+k,), which shows that calculations of current 
using k=0.8 may be made with an error less than that of k only up to 
ro/r=2.55, and above that point the error is greater than that of k, 
becoming infinitely great as 7o/r approaches infinity. 

In the case of spheres the value k=0.8 yields an error of 5.6 per cent 
for r/ro=2.0 where the emitter is inside, and 7.9 per cent for ro/r=2.0 
where the emitter is outside. Thus when current is flowing in either 
direction between any curved surfaces where the radius of the outer 
surface is not more than twice the radius of the inner, the amount of 
current flowing per unit area of a surface four-fifths of the distance from 
the emitter to the collector is very nearly the same as the current between 
parallel planes. In the case of coaxial cylinders with the emitter inside 
this same rule applies with similar accuracy even when the radius of one 
is twenty times the radius of the other. 

The writers wish to acknowledge the assistance given by H. M. Mott- 
Smith in the mathematical analysis of this problem. 
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THE MAGNETOSTRICTION OF A MAGNETITE CRYSTAL 
By C. W. HEaps 


ABSTRACT 


Magnetostriction of an octagonal magnetite crystal along three axes 
was measured in fields up to 5000 gauss. The curves obtained are all similar 
in form, first convex towards the H axis, then concave, tending toward a 
maximum near 5000 gauss. With a longitudinal field of 5000 gauss, the digonal 
axis expands 30(10)~‘ per cent, the trigonal axis 12(10)~* per cent, and the tet- 
ragonal axis contracts about 4(10)~* per cent. With an equal transverse field 
the percentage changes are a contraction of 44(10)~‘, a contraction of 28(10)~* 
and an expansion of 4(10)~‘, respectively. Tests indicated that the effect per- 
pendicular to a plane is independent of the direction of the field in that plane. 
No departure from cubic symmetry was established. To explain these results 
a model of the magnetic element in magnetite is suggested, similar to Ewing's 
recent model. The experiments favor the view that the rotating magnetic 
element is a group of non-parallel electron orbits inside the atom. 

Permeability of a magnetite crystal was roughly determined to vary 
from 1.2 to 1.38 as the field strength increased to 5000 gauss. 


N a recent paper! the writer described some magnetostrictive effects 

obtained with a crystal of magnetite ground into the form of a sphere. 
The experiments, however, were not of a definite enough nature to permit 
the exact correlation of the magnetostriction with the directions of the 
crystallographic axes. The great differences in the character of the 
phenomena for different directions in the sphere made the location of 
these axes with respect to the magnetostriction a matter of some im- 
portance. Further experiments have now been performed and the 
magnetostriction definitely correlated with the axes. 

The apparatus for measuring dimension changes, except for a few 
modifications made in order to facilitate the adjusting of the specimen, 
was the same as described in the previous paper. An electromagnet 
smaller than the one of the previous experiments was used in this later 
work. It had pole-pieces 5.8 cm in diameter and these were kept 3.3 cm 
apart throughout the work. The field, therefore, was probably not as 
uniform as in the case of the larger magnet. Tests with a flux-meter 
and search coil showed the field midway between the pole-pieces to be 
somewhat less than 2 per cent smaller than the field very near the pole- 
pieces. 

The magnetite crystal was an octahedron from Henry County, Vir- 
ginia. It was quite perfect, with principal axes about 1.8 cm long. Before 


1 Heaps, Phys. Rev. 22, 486, 1923 
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making measurements on this crystal the corners and edges were rounded 
off by grinding. It was thus easier to clamp the crystal in the measuring 
device; also the specimen was made more nearly spherical and the in- 
tensity of magnetization presumably more uniform. 

Magnetostriction of the specimen was investigated with respect to 
three different crystallographic axes, the tetragonal, the trigonal, and 
the digonal. Magnetite belongs to the cubic system, the eight faces of 
the octahedron being equilateral triangles. The three fundamental axes 
of the cube, the tetragonal axes, join opposite vertices of the octahedron. 
The four trigonal axes, each of which connects opposite vertices of the 
elementary cube, have the direction of the normals to the faces of the 
octahedron. The six digonal axes are each perpendicular to a pair of 
parallel edges of the octahedron. 
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Fig. 1. Longitudinal magnetostriction along different axes. 


The curves of Fig. 1 represent the magnetostriction obtained in the 
direction of the field for different values of the magnetic field H (measured 
in the absence of the crystal). Here dL/L is the change of length of an 
axis divided by the total length of that axis. The longitudinal magneto- 
striction along the digonal, trigonal, and tetragonal axes varies in much * 
the same way with the field as the latter is increased, but the contraction 
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(plotted below the // axis) along the tetragonal axis is smaller than the 
expansions (plotted above the // axis) which occur for the same fields 
along the trigonal and digonal axes. It is to be noted that two scales 
are used, the curve for expansion of the digonal axis being plotted with 
reference to the scale giving the larger values. 

Fig. 2 gives the magnetostriction along three different axes, the field 
being for each curve normal to the axis specified on the curve. Contrac- 
tions are plotted below the H/ axis, expansions above the axis. It appears 
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Fig. 2. Transverse magnetostriction along different axes. 


that for transverse magnetic fields the tetragonal axis expands while the 
trigonal and digonal axes suffer contractions of considerably larger 
magnitudes. 

The curves of Figs. 1 and 2 are all similar in appearance. They are 
also similar to the curves for a magnetite sphere which are given in the 
previous paper. The curves for the sphere, however, are plotted against 
magnetic induction and are extended to higher values of the field so that 
the existence of maxima is more clearly indicated than in Figs. 1 and 2. 
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The field instead of the induction is plotted in the present work because 
of the difficulty of measuring the latter quantity accurately when the 
specimen is not spherical. Approximately correct values were obtained, 
however, by using a fluxmeter connected to five turns of wire wound 
around the specimen. In this way the permeability of the magnetite 
was found to vary from 1.2 to 1.38 for the field values plotted in the 
curves. 

In its magnetic properties magnetite has been found? to show departures 
from cubic symmetry. To investigate this question the transverse mag- 
netostriction for the different tetragonal and trigonal axes was measured. 
It was found that slight departures from accuracy in the setting of the 
crystal caused considerable changes in the magnitude of the magneto- 
striction so that small deviations from cubic symmetry could not be 
detected easily. It seemed possible, however, that one of the tetragonal 
axes expanded more than the others under the influence of the field. 

From the curves of Figs. 1 and 2, and on the assumption of cubic sym- 
metry, the curves of Fig. 3 were constructed. Here the maximum 
magnetic field is supposed to be rotated in one of the planes of symmetry 
of the crystal, the angle @ between the tetragonal axis and the field being 
plotted as abscissas, the corresponding magnetostriction as ordinates. 
From these curves it appears that a longitudinal field rotated in the axial 
plane (100) produces two maxima of extension and two maxima of con- 
traction in a half revolution, the contractions being small compared with 
the expansions. When the field is rotated through 180° in the (011) plane 
only one maximum of expansion and one of contraction are obtained. 

For the transverse field rotated in these same two planes the magneto- 
striction at right angles to the field is taken in the plane of rotation. 
The curves are similar to those for the longitudinal field except that 
there is an interchanging of contraction and expansion. If a crystallo- 
graphic axis expands in a longitudinal magnetic field it contracts by about 
the same amount when turned transverse to the field; or if it contracts 
in the longitudinal field the transverse field will make it expand by about 
the same amount. 

A number of experiments were made to see whether the magneto- 
striction perpendicular to a given plane varied with the direction of the 
field in this plane. The transverse field was found to produce the same 
expansion of a given tetragonal axis (within the limits of error of the 
experiment) when directed along one of the perpendicular tetragonal 
axes as when parallel with one of the perpendicular digonal axes. Ap- 


? Kunz, Bull. Nat. Res. Council, 3, 177, 1922 
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parently, therefore, rotation of the field in the (100) plane does not change 
the magnitude of the magnetostriction normal to this plane. Further 
experiments made it seem probable that the magnetostriction normal to 
the (011) plane and to the (111) plane is unaltered by rotating the field 
in these respective planes. 


DISCUSSION OF RESULTS WITH REFERENCE TO THEORIES 


In the previous paper it was suggested that Ewing’s recent model of 
the magnetic atom could be made to fit in with magnetostriction phe- 
nomena. The present experiments seem to impose the following con- 
ditions on this model in magnetite: 
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Fig. 3. Magnetostriction for different directions in given planes. 


1. In one complete rotation of the Weber element in the (100), (010), 
or (001) planes it passes through four positions of maximum force which 
tend to expand the crystal and four positions of maximum force which 
tend to contract it. This force is presumably exerted by the Weber 
element either on the outer parts of its own molecule or on neighboring 
molecules. In the former case it must be supposed that the molecule is 
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strained and that the straining of the molecules alters the distances 
between their centers. 

2. The Weber element must possess laterally a mechanism for pro- 
ducing the same type of force changes as are produced longitudinally. 

3. The Weber element, or the rotating mechanism, must be confined 
to the inner part of the molecule. Otherwise a magnetic field, by rotating 
the outer parts of the molecules, could produce changes in those prop- 
erties of the crystal, such as cohesion, etc., which depend on the outer 
parts of the molecules, and the crystallogrpahic axes could have their 
directions changed by means of a magnetic field. 

The following modification of Ewing’s model will satisfy the above 
conditions. Let the fixed outer portion of the molecule be represented 
by a cube with twelve magnets fixed respectively to the center of each 
edge, the north pole of each magnet being directed inwards and the axes 
of the magnets being supposed to lie along the digonal axes of the cube. 
The rotating element may consist of a similar cube centered within the 
first, the magnets of the inner cube having their south poles directed 
outwards towards the north poles of the outside cube. An externally 
applied field must exert a directive influence on this cube; also the 
longitudinal magnetostrictive effects are the converse of the transverse 
effects. Hence it is necessary to assume that one pair of magnets, with 
axes on the same line, have their north poles protruding, and one of these 
two north poles is slightly stronger than the other. 

Fig. 4 gives a diagram of the model. Here the inner rotating cube is 
considerably simplified for the sake of clearness. Only the magnets are 
shown, and each of the poles marked 2 is supposed to take the place 
of four south poles making angles of 60° with each other. The poles of 
the Weber element in the position shown are supposed to lie very close 
to the planes of the outer cube in which the north poles are situated. 
These latter poles are represented by circles. In the diagram, if N’ is 
stronger than N the external field may be supposed directed along Ox. 

Magnetostriction can occur as follows: In the position shown each 
pole marked S will attract those north poles of the plane in which it lies. 
The two & poles will exert attractive forces on all the poles of the ex- 
ternal shell, the compressive force being less along the direction Oz than 
along the other two axes. The poles N and N’, lying approximately in 
the two yOz planes of the diagram, will exert repulsive forces on the 
north poles of these respective planes. Hence all the forces named tend 
to produce contraction in the direction of the field, that is, along the 


tetragonal axis of the crystal. Expansion in transverse directions would 
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result if the repulsive effects of the N and N’ poles exceed the other at- 
tracting forces. 

It seems probable that the action of the Weber element on similar 
elements in neighboring molecules should be considered. In this case 
contraction in the direction of the field would result provided the follow- 
ing condition is satisfied along Ox: compressive force of S poles + com- 
pressive force of = poles > expansive force of N and N’ acting on the 
N’ and N poles of neighboring elements. To produce transverse expansion, 
along Oz for example, we must have for this direction: compressive force 
of S poles + compressive force of = poles < repulsive force of N and 
N'+ repulsive force of = poles on the = poles of neighboring elements. 
In Fig. 4 the inner system is shown in a position which would not 


























Zz 


Fig. 4. Model of the magnetic element in magnetite. 


be stable when forces exist between the different systems. For stabil- 
ity under these conditions the Weber element must be rotated slightly 
about the Ox axis. However it is possible in the rotated position to ex- 
plain magnetostriction in the same way that it is explained above. When 
the field is directed along the digonal axis of the cube a large expansion 
in the direction of the field is to be expected and a similar large contraction 
transverse to the field, because in this case each pole of the Weber element 
is directly under a pole of the outer shell. 

If the influence of neighboring atoms is to be considered in the fashion 
outlined above it is obviously of importance to know the crystalline 
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structure of magnetite in detail. W. H. Bragg* has analyzed this mineral 
and fcund it to have a space lattice similar to that of diamond. The 
(100) planes consist of Fe,O,4 planes having the spacing 2.07 A interleaved 
with Fe planes. The (110) planes consist of Fe,O, planes 2.94 A apart 
interleaved with FeO, planes. The (111) planes are Fe planes 4.8 A apart 
interleaved w:th several iron and oxygen planes. On the average it ap- 
pears that the (110) planes have the largest spacing; it is also these 
planes which experience the greatest normal displacements under the in- 
f'uence of parallel or normal fields. It is possible that the greater spacing 
of these planes gives them greater freedom of motion. The weaker forces 
of the Weber elements on each other across the greater interplanar 
distance does not appear to be of importance in the model suggested 
above, since it is assumed practically negligible in this direction. 

This model of the magnetic atom of magnetite is to be considered 
merely as suggestive. When the magnetostriction of a pure iron crystal 
has been investigated it may be possible, knowing the space lattice of 
magnetite, to devise a much more satisfactory model. The one suggested, 
however, agrees with the concentric shell theory of the atom and it meets 
the demand that the rotating element should be only ‘a small part of the 
atom. On modern theories of atomic structure the magnets of this model 
would naturally be replaced by electron orbits. In this case there is 
objection to the model on theoretical grounds, since for stability it ap- 
pears that the orbits should be coplanar. It is entirely possible that a 
Weber element could be devised with parallel magnets and that this 
element could be made to give the magnetostriction found experimentally. 
However, the assumptions regarding the nature of the forces do not ap- 
pear to be simple enough. The approximate equality of the magnitudes of 
the transverse and longitudinal magnetostrictions suggests similar mech- 
anisms for the production of the effects. It seems natural, therefore, 
to suppose that if aligned north poles produce longitudinal expansion 
then south poles aligned in this same direction produce the transverse 
contraction. Thus magnetostriction phenomena in magnetite favor the 
view that the magnetic element is a group of electrons rotating in orbits 
which are not parallel with each other. This group is probably one of the 
inner shells of the atom. 


Tue Rice INstTITUTE, 
Houston, TEXAS. 
February 4, 1924. 


3 W. H. Bragg, Phil. Mag. 30, 305 (1915) 
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ELECTROMAGNETIC INDUCTION IN A HOMOGENEOUS 
SOLID CONDUCTING SPHERE ROTATING ABOUT AN 
AXIS PERPENDICULAR TO A UNIFORM 
ALTERNATING MAGNETIC FIELD 


By ALEXANDER Marcus 


ABSTRACT 


(1) Mathematical theory. The differential equations of the field are 
solved by a method first given by H. Lamb. The induced currents are found 
to flow in concentric shells; to an observer moving with the sphere they would 
be equivalent to two sets of steady currents revolving in the sphere with fre- 
quencies equal respectively to the sum and the difference of the frequencies of 
the impressed alternating field and of rotation of the sphere. The currents may 
be supposed to diffuse from the surface inward in the usual manner; at higher 
frequencies the’ penetration is smaller, resulting in a ‘‘skin effect’ for high fre- 
quencies. Equations are also derived for the resulting magnetic field and for 
the torque on the sphere and for the logarithmic decrement in the case of an 
oscillating conducting sphere. (2) To test the theoretical results, measure- 
ments of the logarithmic decrement of a solid metal ball oscillating about a vertical 
axis in an alternating field produced by a long horizontal solenoid were made, 
and the results for three balls were found to agree with values computed from 
the theoretical expression for the torque within the limit of experimental error, 
which was not greater than a few tenths per cent. 


HE subject of electromagnetic induction in spherical conductors has 
been investigated by many writers and the special cases that have 
been treated thus far may be divided into two classes, dealing either with 
a sphere at rest in an alternating magnetic field or with a sphere rotating 
in a stéady field. Solutions of these cases have been published by 
C. Niven,! H. Lamb,? J. Larmor,* J. J. Thomson,‘ H. Hertz,®° and R. 


Gans.® 


So far as the writer is aware, no one has hitherto published the solution 
of the problem of a sphere rotating in an alternating magnetic field. 
The general problem of rotation in an oblique field breaks up into two 
cases, for a magnetic field that is unsymmetrical to the axis of rotation 
may be resolved into two components, one that is symmetrical to the axis 
and another that is perpendicular to the axis. The alternating sym- 
metrical component produces the same effect, so far as the currents in 


1 C. Niven, Phil. Trans. part II, 1881 

? H. Lamb, Phil. Trans. part II, 1883 

3 J. Larmor, Phil. Mag. Jan. 1884 

* J. J. Thomson, Recent Researches, pp. 546-557 

* H. Hertz, Collected Papers 

* R. Gans, Zeit. f. Math. und Physik, 48, No. 1, 1902 
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the sphere are concerned, whether the sphere be stationary or rotating, 
for the number of lines of magnetic force that thread through any closed 
curve in the sphere remains unchanged by the rotation. We shall there- 
fore proceed directly to the case of rotation in the magnetic field that 
is perpendicular to the axis. 

Let us suppose that the sphere is rotating with constant angular 
velocity w about the Z axis in an alternating field 7 cos pt parallel to 
the X axis. In this case the induced currents and the magnetic field are 
exactly the same as they would be in a stationary sphere placed in an 
alternating field revolving with angular velocity—w, referred to a set 
of axes X, Y, Z fixed relative to the apparatus producing the field. The 
components of the latter are derivable from the potential 7x cos pt. As 
we shall use spherical coordinates throughout this paper we may write 
this expression as //r sin 6 cos ¢ cos pt. As viewed by an observer 
moving with the sphere, the impressed field is given by the gradient of 
Hr sin 6 cos (go+ wt) cos pt, since the only coordinate that is influenced 
by the motion is g and the relation between ¢ and ¢» is given by g¢=¢+ wt 
where gp pertains to the system of axes revolving with the sphere. 


? 











g=Zreot 


Y 


Fig. 1 


Mathematically the problem is now reduced to the case of a stationary 
sphere in an alternating field whose components with respect to the r, @ 
and go coordinates are [H sin @ cos pt cos (go+wt)]|,, [H cos @ cos pt 
cos (yo+ wt)],, and —[H cos pt sin (go+ wt) |, or {3H sin 6 cos [got (w+ p)] 
+3H sin @ cos [go+(w—p)é]} , etc. To find the induced effects we may 
now employ Maxwell’s equations for stationary media. The following 
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symbols will be used in the equations: 1/, the resultant magnetic in- 
tensivity: E, the resultant electric intensity; 7,6,¢9 the polar coordinates 
of a point; c, the velocity of light; A, the conductivity of the metal; a, 
the radius of the sphere; m, the frequency of the currents in the primary 
coil; m, the frequency of rotation of the sphere. M and E are both 
measured in electromagnetic units. 

The equations pertaining to a particle within the sphere are 

curl M=4nrdE; (1) 
curl E=dM/dt; (2) 
civ E=0=divM ; (3) 

vV°?M=4rd\dM/dt; (4) 
V°E=4rd\dE/dt. (5) 

As we may neglect the displacement currents in this case, we get for 
the region outside of the conductor the equations V*AJ=0 and V?E=0. 

The impressed field gives rise to two types of electromagnetic dis- 
turbances, the free or transient vibrations and the forced vibrations. 
If we allow sufficient time for the free disturbance to subside, we need 
concern ourselves only with the forced vibrations. Moreover the problem 
of determining the precise mode in which the transient effects are dis- 
sipated in time and space in a stationary sphere has been completely 
solved by H. Lamb.’ His solution will apply directly to the case of a 
revolving sphere, for as regards the transient effects the frequency of the 
impressed field is of no importance and as we have replaced the motion 
of the sphere by an appropriate variation of the impressed field, angular 
velocity of the sphere also can not influence the transients. 

The part of the problem which has not been solved heretofore is based 
on the solution of the equation of the type V?M=42\dM/dt which will 
also satisfy div M =0 as well as the boundary conditions. We have already 
proved that relative to the rotor, the impressed field has two frequencies 
proportional to p+ and w—p respectively, where p=27n and w=2rm. 
We may therefore assume that M and E vary as e“ where g=(w+))j 
where j= ./—1 and k?=4mdqg. Ec. (4) may now be put into the form:* 

(rM,) Aa(sin 00M ,/d6) 1 a°M, 
- —+- - —~+4— —- ——-+k?M,=0, 
rer? r? sin 000 r?sin?@ = dg? 
with two similar expressions in Mg and Mg, and div M=0 may be ex- 
panded into 


0 


1 :*.. 1 dMe 
(IT) —— ay — sin 6 Me+— —=0 


r® Or rsin@ 06 rsin@ dg 


7 Lamb, Phil. Trans., part II, 1883 
8H. Weber, Die Part. Diff. Gleich., pp. 336-337. 
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Lamb has long ago given two particular solutions of the systems (1) 
and (II) but only one of these solutions is appropriate to our problem. 
Bromwich® has recently proved that the general solution of Maxwell’s 
equations can be derived by superposing two special solutions which can 
be obtained by taking first E,=0 and secondly 1J,=0. These two solu- 
tions are of a reciprocal character and when EF, =0 none of the components 
of M can vanish, while if 17,=0 none of the components of E can vanish. 
It is obvious that the particular solution which makes M,=0 is inap- 
plicable to this case. 
The appropriate form of the solution is 


2 as 
M,=—}) Ofilkar)—e™; 


Mi= Sf afet eur) — Hf bar) = St eo 
Mi = fe“ T ee > 
* tase , r ao. 


Me= z 2fo( kur) — Rr*fa( bar) | —— OB ant ; 


}rsin @ Ago 


(7) 





2 OSn 
Ee= } nf Safu(kar) =e 


} 

S,=A,r sin 6 e*, where A is a constant. fo(Rar), fi(Rar), and fo(Rar) 
all satisfy 

@f, 2(n+1) dfr 

dx? x dx 


+f.=0, 


where 

x=kar, kk=3V2 (14+j) V4ad(wtp), and ke=}3-V2 (1+7) V4irX(o—p) 

and 
n=j(wtp) and q.=j(w—p). 

Eqs. (6) and (7) were obtained from the corresponding expressions in 
Cartesian coordinates given by Lamb in his third edition of ‘‘Hydro- 
dynamics” on pages 572-573. The differential equation in f, corresponds 
to Lamb’s differential equation in R, on page 478 of the Hydrodynamics. 
It is necessary to use the following recurrence relations between fo, f; and 
fe in order to make (6) and (7) satisfy the field equations 


1 df,(x) dfn 
fnsi(x) = —- ae, and x Dns etl) fam fas. 
d dx 


x x 
* Bromwich, Phil. Mag. (July 1919) 
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Since V?/=0 outside of the sphere, the magnetic intensity due to the 
currents is derivable from a potential. We may therefore write for this 
part of the field 








2 0 " 
ur =-S ot (2) 


1 rsin@ddg\r' ] 







where x=B,rsin 6 é**. 


In order to find the A,’s and B,’s in (6), (7), and (8) we apply the 
boundary conditions for 4. The tangential components of M on both 
sides of the surface of the sphere are equal, and since the permeability 
of the metal is 1, we may also equate the radial components of M on both 
sides of the boundary. These two conditions give us a pair of simultaneous 
equations in the A,’s and B,’s corresponding to each of the two fre- 
quencies. 

At r=a 











—6fi(kia)A, sin 0b Pt" = 2B, sin 6 MPT 4 LT sin 9 AM Meter 
and 














|2fo( kia) — k2a? 2(kia)|Aicos0 fet or 
= —B, cos 0 ef? FPO 4 LT cos 6 i Pt or 

with two similar equations connecting A» with B,. The terms in H on 
the right represent the impressed field. Therefore 

(a) 6fi(kia)Ai:+2B,= —3H; 

(b) [2 fo(kia) —k°a*fo(kia)]A1+B, =3H; 
but 3f\(kia) = —fo(kia) —k,?a*fe(kia), and substituting for k,?a%fe(k,a) in 
(b) we get 

(c) 3[fo(kia) +fi(kia)]A1+B, = 31. 

From (a) and (b) we find 








H —HA\3fi(k 
i se [ fi ( 12) + fo( kia)] 
4fo (kia) 4 fo( kia) 
Similarly we find . 
H —H\3 
A» and B.= fa (haa) + folltea)| 











~ Afo(Rea)” Afo( kea) 
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The value of M due to the currents and the value of E are now given by 
the real components of the following expressions for space inside: 


M,= 3a sin fae 17) ‘iletotod Ft (ker) allot wp | 
fo(kia) fo( koa) 


fo(Rir) +fil kr) cilertiotanl 4 JoChar) + filher) | foslen “a 
fo( kia) fol Rea) 

Mow oa pene cileetrirtond 4 Jolher) + Ail ker) Jletw- ml 
4 fo( kia) fo( koa) 





3 
Me=—Hcos@ \ 
4 














Eg= -j-r cos 0} (wtp yithar) ) 


clototwnl + (yp) SiRer) eilet(w-pyl 


\ 
fo( kia) fol Rea) f 


fi(kir) cilet(otonl pte) ciletw- pul | 
fo(kia) —_ Yea)" f. } 


fo(x) may have any of the forms sin x/x, cos x/x, &*/x, and fi(x) may have 
any of the forms 


3 (10) 
= “Hr} (w+? —— 


1 d sinx 1 d cosx 1d ¢» 


xdx x ' xdx x ' xdx x 
The only form of f:(x) which is applicable to space inside, where x may 
become 0, is the first, and at x=0, fi(x) =1/3. 
The two cases of interest are first the one in which p and w are small 
and ka small, and second, the one in which p and w are large and ka large. 
For small values of x, fo(x)=1 (approx.), fi(x)=1/3. 
For large values of x, fo(x) = &*/2jx, fix) = &*/2x?, 
Substituting in the general formulas (9) and (10), we get for kr small, 
inside: 
Eg=13Hr{ (wtp) cos [go+(w+p)t]+(w—p) cos [eo+(w—p)t}} (11) 
etc. 
M,=—3H sin @{cos [go+(w+p)i]+ cos [yo+(w—p)e]} (12) 
etc. For large values of kr inside, we must take the real parts of the 
following expressions: 


rn. “| ft. —mi(a—1)—j| mia—1) +0 /4] elletoted 


8/7 r 
4 4/2. e— (4-1) — jl mia—) +0 /4) efovwe-ot} (13) 
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i a ee 

M, = ———~ —sin 05 —_—_———- cos| + (w+) 8] 
4Vrr \ VX(wt p) 

e7 m(a—1)—jmila—s) 


\ 
+————_——— cos |got+(w—p)t (14) 
VX(w—p) ae ly 


m=V2rMwtp), m= V2rMw—p) or V27A(p—w) according as p<w, 
or p>w. 

From the expressions for M/ and E it follows that an observer moving 
with the sphere will detect ‘a: pair of revolving fields and two sets of re- 
volving currents. 

Let us now consider the penetration of the disturbance within the solid. 
For small values of and w the current density on any shell varies di- 
rectly as the radius of the shell, being a maximum at the surface and zero 
at the center. The impressed magnetic field being uniform, the same 
inductive effect is produced on all concentric shells. 

For p and w large, the penetration within the conductor takes place 
as a propagation of two highly damped waves inward from the surface. 
At the surface, the current density increases with the conductivity of 
the material, for the components of the current density \E vary as 
V\(w+p). At any depth a—r beneath the surface, the current density 
decreases with an increase in \ owing to the very rapid decrease in 
e~V2*\@+?), The surface of a good conductor therefore screens the in- 
terior against rapidly varying external electromagnetic disturbances 
irrespective of whether the induced effects be due to the rotation 
of the conductor or the variation of the field or to the two com- 
bined. In the case of a conductor rotating in a field that is not 
symmetrical to the axis of rotation the surface can shield the interior 
only against non-symmetrical components, for the rotation in the sym- 
metrical field is not accompanied by induced currents. These facts lead 
to the view that the direct induction creates only superficial currents 
at the start. These currents, were there no resistance, would continue 
to flow undiminished and would prevent the interior of the solid from 
being affected. The sphere may be regarded as made up of an infinite 
number of conducting shells. As the currents of the outer shell decay by 
resistance this shell ceases to be a complete screen. Currents are then 
created in the inner shells by direct induction exactly as in the outer one, 
and thus the whole sphere becomes pervaded. 





THE MECHANICAL REACTION BETWEEN THE CURRENTS IN THE 
SPHERE AND THE IMPRESSED MAGNETIC FIELD 


If F be the mechanical force per unit volume, then (Maxwell’s Treatise, 
vol. II, §603), F=[cxM], the right hand member being the vector 
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product of the impressed magnetic field alone (since the rotor does not 
produce a torque on itself), and the current at the point in question. 
Let X, Y, Z, a, 6, c, u, v, w be the components of F, M and C along the 
Xo, Yo and z axes in the moving system. Then 
ij k 
=|4 v w 
i b ¢ 
and the couple per unit volume is equal to 
i j k 
IrXF]=|xo yo 2 
aA Yaz 
where r=xot-+yoj+sk. The couple around the axis of z is {ff (xp Y—yoX) 
dxdydz; but X=bw-—cv and Y=cu—aw. Substituting we get, 
yoX — xo VY = byw —cyov — cxou + axow = (axot+byot+cz) w—c(uxo+vyo+wsz). 
Since E,=0, wz+vyot+uxo=0, also axo Xbyp+cz=rM, and w=dEz sin 6; 
hence the torque around the z axis is {{fXE,M,r sin @6dS where dS= 
r? sin 6 d@ dg dr or 


r=a O=% e=2n 
Torque ff f f Ee M, rsin?6d@dgdr.+* (15) 
r=0 6=0 ¢=0 


Inserting for /, only the external field, we find for the average resultant 
torque 
af, (kir) filker) 
G = P| (w+) rdrt(o-p) [= nar| (16) 
( 0 fo( kia) fo (Rea) 

For p and w small, 

fi(kir) _Silker) =7 

fo( kia) ~ fo( koa) 3’ 


hence the resultant torque is 
= ~7O Hw. (17) 





In so far as the reaction upon the impressed field is concerned, the cur- 
rents in the sphere flow as though they were confined to tubes. coinciding 
with great circles through the axis of rotation, for in the expression for 
the torque the only component of the current that appears is Eg. More- 
over if we calculate the line integral of Ey about a line whose plane is 
perpendicular to the axis of rotation we find that it vanishes, i.e. 


2e 
fe sin 6dgp=0. 
0 


This must necessarily be true since there is no impressed field parallel 
to the axis of rotation. 
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EXPERIMENTAL CONFIRMATION 


The foregoing theoretical development was tested experimentally in 
the following manner. A metal ball was suspended by means of a torsion 
fibre at the center of a long solenoid carrying a steady alternating 
current. The suspension fibre was perpendicular to the direction of the 
field. The ball was given a start and allowed to oscillate freely in air, 
and an observation of the decrement of the oscillation was taken. The 
field was then thrown on and the decrement was again determined. The 
difference between the two decrements gave the component due to the 
induced currents alone. The latter could be predetermined from theory 
as follows. 

The equation of motion of the ball is given by 


ap dp 
Pe 18 
wea” U8) 


where P is the moment of inertia and k the coefficient of damping includ- 
ing air damping as well as electromagnetic damping. 7 is the restoring 
moment for unit angular deflection of the elastic suspension fibre. The 
solutiop of (18) is given by 

g= ge"? cos (2rt/T) (19) 
where 7 is the period. 

If 5 be the logarithmic decrement of the oscillation we can prove by 
means of (19) that 6=k7/4P. 

In the experiments the period of oscillation of the ball ranged from 
about four to fifteen seconds; therefore it was permissible to compute 
the theoretical value of 6 on the assumption that the angular velocity w 
was constant at any time. 

The first experimental test was performed to check the result given by 
(17) for the average resultant torque. This expression for <Z should be 
strictly true, according to the theory, for impressed fields of zero fre- 
quency. Accordingly, metal balls having different resistivities were sus- 
pended in turn in uniform constant magnetic fields and allowed to oscil- 
late. The damping due to the induced currents was measured in each 
case and the logarithmic decrement for each ball was compared with 
the value given by the theoretical expression. The latter is found in the 
following manner. 6=k7/4P where k=-Zw and, for a sphere, P = 3a2M; 
hence 6= TH*xa*/12Mo, where g is the resistivity of the material or 1/A. 

In every case the values of 6 (the decrement) as given by theory and 
experiment differed by not more than one tenth of one per cent. This 
method was later employed to determine o for copper and for lead from 
the values of 6 obtained experimentally. The values of ¢ obtained from 
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the formula agreed with the published values as closely as the different 
published values agreed among themselves. 

The test of the theory was then extended to those cases for which the 
expression for 6 could be put into the form of a rapidly converging power 
series in the variable »=22aV/f/o where a is the radius of the ball and f 
the impressed frequency. The reason for thus restricting the test was that 
the separation of the real and the imaginary parts of the integral of (16) 
becomes extremely laborious except when we restrict the investigation 
to cases for which v<z. In the latter cases the expansion is made as 
follows: 


—mneomagiGy a 
0 fol Rna) x? 


@ fi(Rar) ace cot x =] 
x 


where x=,a: but 


On substituting in (16) and changing the variable x into v by the rela- 
tion v=2ravVf/o we get 


T = Ps (H*ra*/o)(1—0.11430'—7 .73 XK 10-408... )w, 


6= (TH? ra®/12Mo)(1—0.1143v'—7 .73 X10-5v8...), 


APPARATUS 


The ball was suspended with its center in the middle of a solenoid 
122.4 cm long having an average diameter of 7.1 cm. The solenoid con- 
sisted of two separate but identical parts each wound with four layers 
of No. 16 cotton-covered copper wire. The ball was suspended by means 
of a phosphor-bronze ribbon within a box mounted on leveling screws. 
When the solenoids were placed in contact over the ball the two end 
turns touched at nearly all points. There was a very narrow short slit 
at the top to allow the suspension to pass through. Readings of suc- 
cessive amplitudes were taken by observing the deflection of a beam of 
light on a scale about a meter away. 

Balls of different non-magnetic metals were used. The frequencies 
of the impressed field ranged from zero to 580 cycles per second. In every 
case it was necessary to determine the resistivity of the sample of metal 
used. This was done by keeping the field constant and measuring the 
decrement. The resistivity was computed from 69=TH*ra*/12Mc. 
This expression for 49 follows directly from our theoretical expression for 
do in an alternating field by putting f=0. 
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TABLE I[ 


Comparison of experimental and theoretical values of the decrement p. 








Substance a M T z (effective) Fg p (exp.) p (theor.) 





Copper 0.96cm 33.76g 14.9sec. 1.66 amp. 64.7 1.135 1.133 
Copper 0.53 5.86 4.62 4.00 580 1.013 1.013 
Lead 1.17 79.86 5.24 3.17 62 1.003 1.003 








The effective value of H in the expression for 6 is given by H =30.308i7 where 7 is in 
(effective) amperes. 

In each case the experimental value of p is the mean of a large number 
of observations and the probable error of the mean was in no case greater 
than one part in a thousand. The maximum deviation of a single observa- 
tion from the mean was about five parts in a thousand. 


—S—= 


v 

















O 
— 


Fig. 2. Oscillating sphere suspended inside a solenoid, 122 cm long. 





A comparison of the values of p as found experimentally with those 
computed from the theory shows that within the limits of experimental 
error, theory and experiment give the same values. 

In closing the writer wishes to acknowledge the fact that Prof. M. I. 
Pupin suggested the investigation and that he is particularly indebted 
to Prof. A. P. Wills for a number of helpful criticisms. 
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RATE OF EVAPORATION OF LIQUIDS 


THE RATE OF EVAPORATION OF LIQUIDS 
IN A CURRENT OF AIR! 


By T. B. HINE 


ABSTRACT 


The experiments were made in a large ventilating tunnel at Edgewood 
Arsenal near the end of a 280 ft straight section, 6 ft square. The liquid being 
studied was pumped continuously into a copper pan mounted flush with the 
floor, over the edge of which it overflowed into a lower pan. The amount of 
evaporation in several hours was determined by weighing the pans. Results 
were obtained for chlorbenzene, M-xylene, nitrobenzene and toluene, for winds from 
0 to 15 miles per hr, 49 runs in all. The rate of evaporation is! found to be a 
linear function of the wind velocity, and taking the values of vapor pressure 
given in the literature corresponding to the temperatures of the liquid surfaces, 
the relation suggested by DeHeen in 1891 is found to be correct; the number of 
gram mols evaporated per hour from a liquid surface 60 cm in diam. is given 
by M=(0.1+0.10W)P, where W is the wind velocity in miles per hour and P 
is the vapor pressure in millimeters of mercury. Due to variations in the wind 
velocity and uncertainty as to the exact temperature of the surface in some 
cases, the mean variation of observed values of M from the computed values 
is about 10 per cent, but there is no evidence of a systematic deviation of the 
results for any one liquid, although the vapor pressures very from 0.3 to 30 mm. 


I. INTRODUCTION 


HE persistence of a chemical warfare agent in the field is mainly 

determined by its rate of evaporation under the conditions of dis- 
tribution employed. If we know the persistence of one substance for a 
given mode of distribution, and the relation between the physical prop- 
erties of a substance and its rate of evaporation, we should be able to 
predict the persistence of a new substance, and perhaps gain an approxi- 
mate idea of the concentrations in the air over an evaporating area. For 
this reason, the following study of the laws of evaporation was under- 
taken. 

II. HisTORICAL AND THEORETICAL 


The importance of water evaporation phenomena for the agriculturist, 
engineer, and meteorologist has resulted in many investigations and the 
collection of a considerable mass of data, but unfortunately the greater 
bulk of the literature is concerned with the statistics of evaporation, and 
is not of immediate value in determining the physical laws of evaporation. 
Another difficulty that we encounter in analyzing much of the data of 
water evaporation is that the rate of evaporation is very dependent upon 
the humidity of the air which is seldom measured very accurately and 


1 A contribution from the Chemical Warfare Service, War Department. 





80 T. B. HINE 


which is subject to rather rapid and irregular changes, thus introducing 
an uncertainty that necessitates a very cautious use of the conclusions 
based on water evaporation in our study of the rate of evaporation of 
other liquids. 

There is no comprehensive and critical treatise on this subject, but 
fortunately, the several hundred widely scattered publications up to 
the year 1909 have been covered in an excellent annotated bibliography 
by Mrs. Grace L. Livingston,* and the literature since that date may be 
readily found in the abstract journals. An inspection of this voluminous 
literature reveals the rather surprising fact that almost no work has been 
done with liquids other than water and that very few experiments have 
been sufficiently controlled to indicate the effect of each of the factors 
that together determine the rate of evaporation of a liquid in a current 
of air. 

There have been a number of empirical formulas proposed to express 
the rate of evaporation of water, but there is none that has proved suf- 
ficiently satisfactory to be generally adopted, and they have not been 
tested in their application to liquids other than water. 

Since John Dalton’ nearly all investigators have agreed that the rate 
of evaporation of water is proportional to the saturation deficiency. 
This may be expressed as the difference between the vapor pressure of 
the water at the temperature of the water and the partial pressure of 
water vapor in the air. The partial pressure of water vapor in the air 
(or the absolute humidity) is generally determined at a great distance 
from the surface, or before the air passes over the water surface whose 
evaporation is being considered. Apparently assuming that the tem- 
perature of an evaporating water surface would be the same as that of a 
wet bulb thermometer, certain writers following Tate* have expressed 
the rate of evaporation as proportional to the depression of the wet bulb 
below the dry bulb thermometer. With the above assumption this would 
follow from the ordinary psychometric relation that the difference be- 
tween the tension of saturated vapor at the temperature of the wet bulb 
and the actual tension of vapor present in the air at the temperature of 
the dry bulb is proportional to the depression of the wet bulb reading 
below the dry bulb reading. Aside from the undesirability of introducing 
into a physical relation factors dependent upon apparatus construction 
and operation, the main objection to setting the rate of evaporation 


? Mrs. Grace L. Livingston, Monthly Weather Review, June, September, November, 
1908, and February, March, April, May and June, 1909 

3 John Dalton, Mem. Lit. Phil. Soc. 5, 535-602 (1802) 

* Thomas Tate, Phil. Mag. 23, pp. 126, 283, 494 (1862); 25, 331 (1863) 
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proportional to the psychrometric difference is that it takes no account 
of the temperature of the evaporating liquid which may be very different 
from that of a ventilated wet bulb. 

In considering the rate of evaporation of liquids other than water the 
air in most cases will not contain any of the vapor of the evaporating 
liquid before passing over its surface or at a considerable distance above 
the surface, and consequently we shall not be troubled by the fluctuations 
in humidity which so profoundly influence the rate of evaporation of 
water in the open. For any given liquid evaporating into air previously 
devoid of its vapor we might expect that, with the other factors constant, 
the weight evaporated per unit time would be proportional to the vapor 
pressure at the temperature of the liquid. This is equivalent to consider- 
ing that the same volume of air saturated with vapor is removed per unit 
time with other factors constant, since the weight of a substance in a 
given volume of saturated air is, of course, proportional to the vapor 
pressure of the substance. If we now wish to compare the rate of evapora- 
tion of different liquids under otherwise constant conditions, we might 
expect the same volume of saturated air to be removed per unit time in 
each case, and, of course, the weight of substance contained in this volume 
in each case is proportional to the vapor pressure and the molecular 
weight of the substance. Asa result of experiments on the rate of evapora- 
tion of liquids from a bottle in a current of air, DeHeen® concluded that 
the weight evaporated per unit time is proportional to the product of 
the molecular weight by the vapor pressure. This may be expressed in 
another way by saying that the number of mols of liquid evaporated per 
unit time is proportional to its vapor pressure. 

It was for the purpose of testing this relationship in the case of a free 
liquid surface in a horizontal current of air that the experimental work 
described in the rest of this report was undertaken. We may anticipate 
the results of this experimental work, and state that within the experi- 
mental error we found the relationship correct; the number of mols of 
a liquid evaporated per unit time from a given surface is proportional 
to the vapor pressure of the liquid and also to the wind velocity. 


III. EXPERIMENTAL‘ 
VENTILATING TUNNEL AND APPARATUS USED IN EVAPORATING 
EXPERIMENTS 
In order to have the wind velocity under control, the evaporation ex- 
periments were all conducted in one of the large ventilating tunnels 
5 P. DeHeen, Bul. Acad. Sci. de Belgique, 21, 11-25, 214-9, 798-810 (1891) 


® Most of the experimental work here described was performed by Mr. D. Clark with 
general supervision by the author. 
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available at Edgewood Arsenal. The tunnel proper is about 600 feet long 
and parallel with the ground level. It is constructed of wood double- 
walled, with a layer of tar paper between, and painted with asphaltum 
paint. It contains four right-angle bends and is ventilated by four fans, 
driven by four 20 hp a.c. motors, single speed, which give an air velocity 
varying from one mile or less, by the use of swinging gates, up to 40 miles 
per hour, with all fans going and all the control gates open. A low wind 
velocity can readily be obtained that is constant to within about 1/2 
mile per hour. More constant velocities cannot be obtained because 
outside weather conditions are a controlling factor. The apparatus was 
set up in a straight section of the tunnel about 280 feet in length with 
cross section 6 feet square. The wind had a straight sweep of 250 feet 
before striking the exposed pan containing the evaporating liquid. 

The liquids were evaporated in a circular copper pan, 59.06 cm in 
diameter by 1 cm deep, and were pumped continuously into this pan 
from which they overflowed the edge into a larger pan containing the 
bulk of the evaporating liquid. A small pump kept the liquid circulated 
from the lower to the upper pan at a speed which kept the upper pan 
overflowing continuously, insuring a full surface of liquid to the wind at 
all times. The lower pan receiving the overflow was larger in diameter 
by 5 cm and was not exposed to the effect of the wind. A circular hole 
27 inches in diameter, was cut in the floor of the tunnel. A steel plate, 
30 inches square, 1/4 inch thick, containing a circular hole cut slightly 
larger than the diameter of the upper pan but smaller than the diameter 
of the lower one, was placed over the hole and fastened to the floor. 
This brought the edge of the upper pan flush with the surface of the steel 
plate when in position for evaporating. The edge of the lower pan rested 
firmly against a double gasket of hard and soft rubber on the under side 
of the steel plate. This insured an air-tight fit between the lower pan and 
the steel plate. 

The small circulating pump was made of aluminum, eccentric type 
with brass blades, which delivered about 10 cc per minute. The shaft 
speed was 96 r.p.m. driven from a counter-shaft by a small electric 
motor, belt connected. The liquid was discharged just above the evaporat- 
ing surface; it could then be noticed if the pump was functioning properly. 
The driving apparatus, including the under side of the pan, was housed 
directly under the tunnel in a suitable manner. The pan was lowered 
and raised to the tunnel floor by a rope connected to the center of the 
pan by a hook. The other end of the rope extended over two pulleys 
and down under the tunnel to a counterweight. This arrangement was 
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found to be most convenient in lowering and raising the pan without 
spilling the contents. 

A Troemner solution balance No. 80 was used for all weighings of the 
pan after run No. 12. Before this, several types of balances were tried 
but the Troemner solution balance No. 80 was found to be the most 
satisfactory. It was of 20 kg capacity and, when tested for sensitiveness 
at full load of 8 kg, responded to two grams in each case. In our calcula- 
tions we have assumed +4 grams to be the nearest that the beam could 
be read in actual work. 


THE Liguips UsED IN EVAPORATION AND THEIR VAPOR PRESSURES 


The liquids available in sufficient quantities for these experiments 
were as follows: toluene, chlorbenzene, m-zylene, nitro-benzene. All were 
of good quality, but they were purified by us before being used for 
evaporation. 

Toluene (mol. wt.=92.1) was freed from water with fused calcium 
chloride and used without further purification. Its boiling point was 
constant at 110° C and it was considered high grade material pure enough 
for our purpose. 

The vapor pressure at the temperature of the evaporation experiments 
was calculated from the experimental data of Kahlbaum,’ Barker,’ and 
F. Neubeck® by the Kirchoff-Rankine-Dupre formula, log P= —A/T 
—B log T+C, where P is the vapor pressure, 7 the absolute temperature 
and A, B, and C are constants. 

The constants of the equation were evaluated from the pressures given 
by Kahlbaum for temperatures 298.8°, 304.6° and 318.2° K. The equa- 
tion used was: log P= —2075/T—0.9286 log 7+10.7208. Vapor pres- 
sures calculated from this equation agreed satisfactorily with the meas- 
ured values of Barker and Neubeck, for temperatures 273°, 287.8°, 334.2° 
and 336.2° K. 

Chlorbenzene (mol. wt. = 112.53). This was freed from water with fused 
calcium chloride and subjected to fractional distillation through a Hempel 
column containing glass beads, to determine its boiling range only. From 
a total of 500 cc, 465 cc distilled over between 131° and 132° C and 35 cc 
came over between 132 and 134° C using a calibrated Anschutz ther- 
mometer. This means 93 per cent distilling between the range of 
131°-132°C. The boiling point is given in the literature’® as 131.7°C at 


7 Kahlbaum, Zeit. Phys. Chemie, 26, 603 (1898) 

8 Barker, Zeit. Phys. Chemie, 71, 235 (1910) 

®* F. Neubeck, Zeit. Phys. Chemie, 1, 656 (1887) 

10 Ramsay and Young, Jour. Chem. Soc. (London) 47, 642, (1885) 
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756 mm pressure. It was considered unnecessary to distil the amount 
used for the evaporation, so it was used after being freed from water. 

The vapor pressure values were taken from the work of S. Young." 
These were considered reliable, and were interpolated by plotting the 
reciprocal of the absolute temperature against the logarithm of the 
observed vapor pressure and drawing a straight line through the plotted 
points. 


M-xylene (mol. wt.=106.12). This liquid always contains a slight 


amount of the ortho and para isomers which cannot be separated by 
fractional distillation on account of their close boiling points. The ma- 
terial on hand, however, was quite pure and when fractionated through 
a Hempel column distilled between 139° and 140° C. The total amount 
was fractionated before being used for the evaporation experiments. 
It was the only liquid available at the time having a vapor pressure 
between that of chlorbenzene and nitrobenzene. The only vapor pressure 
values found in the literature were those determined over a short range 
by F. Neubeck.* These data were plotted and extrapolated to the lower 
temperatures on a straight line curve. 

Nitrobenzene (mol. wt.=123.08). This liquid was quite pure but was 
slightly brown in color, so it was freed from water with calcium chloride 
and distilled through a straight-glass column. It had a. constant boiling 
point of 210° C throughout and its analyses showed a purity of 96.06 
per cent. The boiling point given by W. Louguinine'? is 210.6° C at 
760 mm. Since the vapor pressures of nitrobenzene as determined by 
Kahlbaum,’ Neubeck,? and Marcellin" do not extend to room tem- 
perature and do not agree in the lower range measured, P. D. Watson 
of this laboratory measured the vapor pressure very carefully by the 
boiling point method from atmospheric pressure down to 25 mm and 
obtained an excellent check with the data of Kahlbaum at high pressures 
although deviating slightly at the lower pressures. This check is of 
considerable interest to us since we have used his values for several of 
the other liquids employed in the evaporation experiments. 

The vapor pressures at the temperatures of the evaporation experi- 
ments were calculated from the equation log P = —3177.307/T — 3.62491 
log 7+19.18057 in which the constants were evaluated from the follow- 
ing data of Watson for the vapor pressure at 210.7°, 197.7° and 103.0° C: 
759.6, 552.0, and 24.8 mm, respectively. This equation reproduces the 
other data of Watson and that of Kahlbaum above 25 mm. Some ap- 

1S. Young, Jour. Chem. Soc. (London) 55, 502 (1889) 


2 W. Louguinine, Ann. Chim.-Phys. (7) 27, 119 (1902) 
8 Marcellin, Comptes Rendus 158, 1420 (1914) 
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proximate measurements of the vapor pressure at room temperature were 


made by Watson-using the method of Menzies." These measurements 
give values that are too high below the room temperature (27.8° C) due 
to a slight amount of nitrobenzene in the connecting tubing out of the 
bath, but the value 0.525 mm at 30.1° C is probably fairly accurate and 
it is seen to agree very well with the value 0.504 mm calculated from the 
equation above. Press of other work prevented more accurate measure- 
ments at the time but these indicate that the values for the vapor 
pressures at room temperature calculated by the equation based on meas- 
urements at higher temperatures are probably very close to the true values 
although their exact accuracy cannot be stated. 


METHOD OF CONDUCTING MEASUREMENTS 


With the evaporating pan in position, the suction fans were started 
and regulated to the desired velocity requiring about 20 minutes to at- 
tain constant readings on the anemometer. The circulating pump was 
started, the thermometer set in place, and part of the weighed liquid 
poured into the pan in such a manner as to cover the whole surface at 
once without splashing. The exact time was noted at the start of evapora- 
tion, and the rest of the liquid then carefully added to avoid spillage. 

The temperatures, anemometer readings and relative humidity read- 
ings were then recorded. The readings were taken as often as was 
thought necessary, depending upon the fluctuations in temperature of 
the liquid, i. e. in some cases, the temperature variation was so great 
that a reading was recorded every 10 or 15 minutes and in others, read- 
ings were required only at 30 minute intervals. The length of time 
for a run depended entirely upon the vapor pressure of the liquid and the 
capacity of the pans. The wind could be switched off at the exact time 
required at the end of each experiment by a remote control switch in- 
stalled under the tunnel. 

When the final temperature of a run had been recorded, a small check 
valve was released in the center of the upper pan, allowing it to drain 
somewhat. This facilitated handling and weighing of the pans, partly 
full of liquid. This was usually done about five minutes before time to 
stop the experiment. The two operators were then in readiness to take 
down the pan at the instant the fans were cut off. The pump was dis- 
connected and the pan lowered by the rope, until the sliding shaft rested 
at its lowest position. The key was pulled out of the shaft, releasing the 
pans, etc. The operators, one on each side of the pan, gently lifted the 


4 Menzies, Jour. Amer. Chem. Soc, 24, 2218 (1920) 
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pan off the shaft and very carefully placed it on the balance. From the 
time the wind velocity was stopped till the actual weight was recorded 
on the balance, an interval of two to three minutes elapsed; therefore a 
slight correction might seem necessary on account of the evaporation 
of the more volatile liquids into still air at this prevailing temperature. 
In runs 49 and 50 an approximate idea of this correction was experi- 
mentally determined for toluene and in run 58 for m-xylene by evapora- 
tion in still air. The greatest rate thus obtained is only about 4 grams 
per minute, hence this loss is within the experimental error of weighing 
and no correction was made for it in the tabulated data. 

This procedure was strictly adhered to in all cases, with the exception 
of some of the first runs where the liquid used was contained in more 
than one bottle, thus making the experimental error greater by 3 or 5 
extra weighings. Another instance was that of adding more liquid during 
evaporation. This addition of new liquid sometimes caused a considerable 
rise in temperature, as noted in the data. 

The total possible weighing error in the weight of liquid evaporated 
is + 16 grams using the method of procedure described above. Assuming 
this error for all runs after No. 12, the per cent error in weighing has been 
calculated for each run and is recorded with the other data in Table I. 

Of the numerous difficulties encountered, the method of weighing the 
pan after evaporation caused the most serious trouble, but as the oper- 
ators became more skillful in the technique of handling the apparatus 
full of liquid, this difficulty was overcome and few losses occurred by 
spillage. 

Throughout the experiments, the temperature of the liquid was kept 
as uniform as possible by circulating the liquid with a small rotary pump. 
At all velocities, a ripple was noticed on the surface of the liquid. The 
main factor during the experiment was that of holding the air to a constant 
velocity. At the rate of 6 miles per hour or less, it was held to within the 
maximum variation of 0.5 mile per hour. At higher velocities, it could 
be held much more constant. It was noticed that outside weather con- 
ditions caused an appreciable effect at the lower rates, primarily due to 
the direction and force of the outside wind at the intake openings of the 
tunnel. The varying speeds of the motors driving the fans also influenced 
this to some extent. 

During the first experimental runs, the upper pan became slightly 
warped and no longer sat perfectly level. The corner of a steel plate 
was placed on the edge of the pan to keep it level during subsequent runs. 
We did not realize at the time that this plate lying close to the pan might 
have some effect on our results. It is possible that the variable turbulence 
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caused by the air flow striking the plate may be the cause of the dis- 
crepancy in some of the data. It is also possible that air leakage may have 
occurred in some cases between the pan and the gasket thus making 
another error to contend with in our measurements. 

Measurement of wind velocity. A Robinson cup anemometer with an 
e'ectrical contact for every 88 feet was used in the earlier work but was 


not of sufficient accuracy for our purpose at the low velocities, although 
it gave very reliable results at velocities above 6 miles per hour. 


For the low velocity measurements, i. e. 6 miles per hour or less, a 
special vane type anemometer made by A. Stroppaine and Company 
of Berne, Switzerland, was used. The vanes are made of aluminum and 
the dial which reads in meters projects on an arm. The instrument was 
calibrated by the makers on a rotating apparatus 10 meters in diameter 
taking the wind of rotation into consideration. A calibration certificate 
accompanied the instrument. This was an extremely sensitive instru- 
ment and all the readings were taken with it on the windward side of the 
pan about 2 inches away from the edge on the steel plate. The axis of 
the vanes was about 2 inches above the surface of the liquid. 

Measurements of temperature. Temperature measurement in the pan 
was taken at the time noted in the experimental data with a 14 inch 
calibrated thermometer No. 7334, reading in 10ths of a degree and having 
a range of —10 to +100° C for 3 inch immersion. It was laid almost 
horizontal across the pan through the hook in the center, the evaporating 
liquid reaching the 3 inch immersion mark at an angle. 

During the first part of the work, a calibrated thermometer was not 
available, therefore an ordinary chemical thermometer range —5 to 
+100° C was used in runs 0 to 14. As it was then broken, we were unable 
to calibrate it and correct the temperature for these runs. An alcohol 
thermometer of the required range was used in runs 14 to 19. This was 
afterwards calibrated against the above instrument No. 7334. The cor- 
rections was applied to the experimental data before they were inserted 
in the tables. 

The temperature of the air in the tunnel was taken by an ordinary 
uncalibrated chemical thermometer, ranging from —5 to +100° C, 
reading in single degrees. 

Measurements of relative humidity. The relative humidity was measured 
with a Tycos hygrodeik and checked at intervals by a sling psychrometer. 
The wicks of the wet bulb were changed once a month and fresh water 
was used every day. The relative humidity was not recorded at regular 
intervals but was taken 2 or 3 times during a run, and in some instances 
5 or 6 times, when the moisture content of the air fluctuated. 
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Barometric pressure. Barometric pressure inside the tunnel was taken 
with a standard Fortins Barometer, U. S. Weather Bureau type, with a 
vernier scale reading to 1/10 mm. The pressures are not included in the 
data, being recorded in only a few runs to determine the pressure change 
for wind velocities of 6 and 12 miles per hour. The effect is small, amount- 
ing to 1/2 mm from normal pressure at 6 miles per hour and an average 
of 1 mm at 12 miles per hour. 


EXPERIMENTAL DATA 


The essential data for 49 evaporation runs have been condensed 
into Table I. Some of the runs were spoiled by accidents, generally 
spillage of the liquid, and the data for these have been omitted. 
The air temperature was generally fairly constant during a run, and it 
has not seemed worth while to record all the individual temperature 
readings, only the average value being given in the table. The liquid 
temperature has such an important influence on the rate of evaporation 
that all the individual corrected temperature readings were used in cal- 
culating the weighted mean vapor pressure, although only the average 
temperature of the liquid is included in Table I. Using the vapor pressure 
data already discussed, the weighted mean vapor pressure has been cal- 
culated for each run as follows: The vapor pressure corresponding to each 
individual liquid temperature reading has been multiplied by the num- 
ber of minutes elapsed until the next temperature reading and the sum 
of these products for each run has been divided by the total time of the 
run in minutes. This weighted mean vapor pressure is recorded in Table I. 

In the runs before No. 14 the accuracy of the weights is not known, 
but it was undoubtedly much less than in the later experiments. The 
thermometer used in this work was not calibrated and the cup anemometer 
was used. These first twelve runs, therefore, should not be given as much 
importance as the later ones. 

As mentioned in the description of the tunnel, the wind velocity could 
only be held constant within about 0.5 miles per hour and although the 
wind velocities recorded in the table are averages of several readings, 
regular readings were not made and, consequently, a weighted mean 
wind velocity cannot be given as in the case of the vapor pressures. There 
is, therefore, a considerable percentage error attached to the wind veloc- 
ities recorded in the data, especially for the lower wind velocities. 


IV. DiscussIon oF DATA 


Since the object of this experimental investigation was to determine 
the relation between rate of evaporation and vapor pressure, the ratio 
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of the mols evaporated per hour to the vapor pressure in mm of mercury 
has been calculated for each run and included in Table I with the data. 
It is at once apparent that this ratio increases with wind velocity and, 
consequently, in order to compare the different runs with each other, 
the ratio has been corrected for wind velocity as follows: In the accom- 
panying figure the ratios of mols evaporated per hour to vapor pressure 
have been plotted against the wind velocity for all the runs. Inspection 
of the figure shows that while the experimental errors cause a considerable 
scattering of the points, the relation between the wind velocity and the 
ratio of rate of evaporation to vapor pressure may be fairly well repre- 
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Fig. 1. Ratio of evaporation to vapor pressure as a function of wind velocity. 


sented by a straight line. The principal experimental error for most of 
the points is due to errors in the wind velocity measurements and this 
is especially true for the low wind velocities. For the runs recorded as 
zero wind velocity the fans in the tunnel were not run but there was 
undoubtedly some air movement and consequently more evaporation 
than would correspond to actuai still air. The line has been drawn to 
give a less value for the ratio at zero wind velocity than the experimental 
points and to represent the mean of the other points as well as may be. 
It would have been possible to represent the mean of the experimental 
values somewhat more closely with a curve but this is hardly justified 
by the accuracy of the data while the simplicity of the straight line 
relationship recommends it for use as a first approximation. 

The straight line in the figure represents the equation: mols/hour 
=(0.1+0.1 wind velocity in miles/hour)(vapor pressure in mm). 
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Using this equation the mols per hour evaporated have been calculated 
for all the experimental runs from the wind velocity and vapor pressure. 
If we omit the runs for zero wind velocity which are unreliable for the 
reasons already discussed, the average per cent deviation for the remain- 
ing 46 runs, neglecting signs, is 10.72 per cent. The calculated rate of 
evaporation is greater than the observed for 17 runs, and less than the 
observed for 29 runs. This agreement between calculated and observed 
rates of evaporation is remarkably close when we consider the errors 
involved in the measurements and that we are concerned with 7-fold 
changes in wind velocity and variations in vapor pressure of 100-fold. 
Such an agreement is adequate proof that the fundamental relation 
underlying the calculation is sound and we may conclude as a first 
approximation that the number of mols of any liquid evaporating 
into a current of air is proportional to the vapor pressure of the liquid. 
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THERMAL DIFFUSIVITY AND CONDUCTIVITY 
OF SOME SOIL MATERIALS 


By L. R. INGERSOLL AND O. A. KoOEpp 


ABSTRACT 


Thermal diffusivities of certain soil materials, quartz sand, sandy clay, 
calcareous earth and packed snow, have been determined by a method based on 
the principle of linear heat flow in a slab whose faces are suddenly chilled. The 
material was packed in a narrow sheet copper box having broad faces, with a 
resistance thermometer in the center. After standing at room temperature for 
two days it was placed in stirred ice-water and the time required for the tem- 
perature to fall half-way to zero noted. This allowed direct calculation of the 
diffusivity, and from a knowledge of the density and specific heat, of the con- 
ductivity. 
HILE the thermal conductivity of a great variety of soils and rocks 
has been measured, there are relatively few diffusivities listed, in 
spite of the fact that in practically all cases of heat flow, save the steady 
state, it is the latter quantity which directly determines the temperature at 
any point.- Moreover, save in relatively few instances,' such diffusivities 
are not the result of direct determination but are calculated from con- 
ductivity measurements by the relation h*=k/cd, where c is specific heat 
and d density. The diffusivity #? in the present case was measured as 
follows. The material was packed into a rectangular box of sheet copper, 
17.5 cm high, 19.5 cm wide and 8 cm thick. In the center was a resistance 
thermometer of No. 38 enameled copper wire wound on a piece of mica 
2 by 4'cm and covered with marine glue, placed parallel to the broad 
faces. This was connected with a Leeds and Northrup recorder. To 
make the heat flow strictly linear and normal to the broad faces, the 
three narrow sides were lined with paraffined wood 2 cm thick, with a 
small air space next the copper. The top had a close-fitting cover. 
After packing, the box was allowed to stand in a pail of water in a room 
at reasonably constant temperature for two or three days; or until the 
recorder showed a long line of steady temperature—aboyt 24° C save 
for snow, which, of course, was initially below zero. It was then plunged 
into a large jar filled with cracked ice and water, and this instant was 
marked as the zero of time on the recorder sheet. In a few minutes the 
temperature at the center would begin to fall and within an hour, in most 
cases, would drop half-way down to its final value of 0° C. The fall was 


1 For instance, Y. Tadokoro, Tohoku Univ. Sci. Reports, 10, 339 (1921) 
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gradual enough, however, so that there was no appreciable lag on the 
part of the recorder. During this time the ice water was vigorously 
stirred, keeping at zero the faces of the box, and therefore of the two 


surfaces of soil. 

It will be seen that we have approximated here the ideal case of linear 
flow in a slab, initially throughout at the same temperature, whose faces 
are suddenly lowered to and held at zero. The temperature at the center 
is given by the equation :? 

@= oo) 1 _4 eW Pry 4 oben ot , 


T 3m 
where /=8 cm in this case. For 0/@)=1/2, h?x*t/i?=0.935 and h? is 
calculated accordingly. Values of h? calculated on the basis of 6/@)=1/3 
and 2/3 were in agreement with those from’ @/0)=1/2. 
The following table lists the results. The conductivities have been 
calculated from the diffusivities. The densities were measured in all 
cases, but the specific heats were as far as possible taken from tables. 








; Percentage Diffusivity Specificheat Density Conductivity 
Material moisture h? c d K 





Quartz sand, me- ’ 
dium fine, dry. . 0.0020 0.19 1.65 0.00063 
Same, with water = 0.0034 0.24 1.75 0.0014 
Sandy clay S. 0.0038 0.33 1.78 0.0022 
Calcareous earth, 
water saturated 3. 0.0019 0.53 1.67 0.0017 
Snow, densely 
0.0041 0.50 0.54 0.0011 








Where comparison of these values with existing data is possible, the 
agreement is satisfactory. 


DEPARTMENT OF PHYSICS, 
UNIVERSITY OF WISCONSIN, 
March 4, 1924. 


? Ingersoll and Zobel, Math. Theory of Heat Conduction, p. 108. 








BOOK REVIEWS 


BOOK REVIEWS 


Foibles and Fallacies of Science, by DANIEL W. HERING—This book will be appre- 
ciated by men of science because of the entertainment which it affords. It can be 
recommended to readers who make no pretense to a knowledge of general science or to 
the mastery of any special field thereof, because of the warnings contained in the book to 
beware of imposters and would-be scientists. Notwithstanding the activities of our 
scientific men, the great masses of our people have little or no comprehension of the 
method of science. Opposition to scientific medical research is assuming alarming pro- 
portions. A large part of our population would favor legislation against the teaching of 
evolution. Rain makers still command a comfortable income. What is the cause of 
this mental flabbiness and this obscurantism in our population? Perhaps the opposition 
of modern educational theory to rigorous scholastic standards, particularly to the serious 
study of mathematics, physics, chemistry and foreign languages in secondary schools, 
may account in part for this deplorable attitude of the popular mind. Any publication 
such as Hering’s book, which contributes to more general enlightenment, is to be wel- 
comed. 

Among the topics treated are astrology (ancient and modern), almanacs making 
weather predictions, transmutation of metals, various perpetual motions including the 
Keely motor, attacks upon Newton's law of gravitation and the wave theory of sound, 
the devining rod, palmistry, Symmes’ theory that the earth is hollow, the moon hoax of 
1835, the Cardiff giant, various prophesies, the actions of charlatans, the elixir of life 
and the fountain of youth. The book is illustrated and contains a short bibliography 
at the end of each chapter. 

Passages in Hering’s book which relate to the history of science are as a rule carefully 
drawn, but a few statements occur which should not be taken seriously. The sentence, 
“at the time of Newton and for more than a century after, the corpuscular theory of 
heat and light prevailed,’’ would be nearer the truth if “at the time of Newton’’ were 
omitted. During the seventeenth century many of the leaders (including Descartes, 
Amontons, Francis Bacon, Newton, Boyle) held that heat was due to motion. As to 
light, Huygens’s and Hooke'’s wave theory had adherents. Newton himself advanced 
the emission theory only with great hesitation, as a second hypothesis, and never defi- 
nitely discarded the wave theory. Newton was able simultaneously to entertain two 
di tinct hypotheses of light. —Pp. XIII+294. Price $2.50. Van Nostrand, 1924. 

FLORIAN CAJORI 


Galvanomagnetic and Thermomagnetic Effects. L. L. CAMpBeLL.—This book is 
another of the series of Monographs of Science edited by Sir. J. J. Thomson and Professor 
Frank Horton. In it is given an historical, experimental and theoretical account of the 
Hall effect and the group of phenomena related to it, including the change of electrical 
resistance, of thermal conductivity and of Peltier effect, in a magnetic field. There is a full 
account of the different methods by which these different phenomena and their variation 
with temperature and with the intensity of the magnetic field have been studied. To this is 
added a complete compilation and a carefully prepared analysis of the data which have 
been obtained on metals, alloys, compounds, electrolytes and gases. Excellent tables 
give in compact form the magnitude of these effects for different temperatures and 
for different strengths of magnetic field. The brief statements of the theories of these 
effects sei ve to show that the existing theories are insufficient and tentative guesses after 
the truth. At the end is a complete bibliography of the literature bearing on the subject. 
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The bringing together of the data and results in this field and arranging them in a 
logical sequence is a valuable contribution, particularly as no detailed account of this 
group of diverse phenomena had previously been published. The importance of this 
compilation in connection with the formulation of any theory of the solid state must be 
recognized. It should be in the hands of all those interested in the electron theory of 
metals.—Pp. XII+311. Price $5.25. Longmans, Green and Co., 1923. A. W. Smita 


Friction, by T. E. Stanton.—The aim of this book, in the words of the writer, is 
“the collection in a single volume of the results of modern investigations into the nature 
and laws of friction.” It is the only book written which is devoted entirely to friction 
and which covers the subject in a comprehensive way; and its author is a recognized 
authority on this subject, well qualified to write such a book. The subject matter is 
very similar to that presented by the author in his article on “Friction” in The 
Dictionary of Applied Physics. 

The sections on fluid“friction, to which more than half of the book is devoted, are 
of particular scientific value. After an outline of the derivation of the equations of 
motion of a viscous fluid a study is made for certain simple cases of the critical condi- 
tion of flow in which turbulence arises. Experimental data are given on turbulent 
motion for certain cases in which the general law of surface friction, which is believed 
to hold, receives an interesting confirmation. The important theories of lubricated 
surfaces are presented with accompanying experimental studies. The remainder 
of the book is taken up with the subjects of solid friction, rolling friction and the trans- 
mission of heat to a fluid passing over a heated solid surface. 

The material, though necessarily rather condensed for the most part, is clearly pre- 
sented. Several small typographical errors in the more mathematical portions require 
correction.—Pp. 180, Longmans, Green and Co., 1924. ARTHUR L. KIMBALL, JR. 


Modern Gypsies. MAry CREHORE BEDELL.—The account of a twelve thousand mile 
motor camping trip encircling the United States, made with his family by Frederick 
Bedell, for many years editor of the Physical Review. Although it does not deal with 
physics, the book will have a personal interest for many readers of this journal. It will 
also appeal to all lovers of out-of-doors, of travel and adventure, for in a simple, straight- 
forward style it gives a continuous series of well selected and interesting details of the 
people, local customs, and landscape encountered.—Pp. 262, 35. illustrations. 
Brentano’s, 1924. 
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CORRECTIONS 


FREE CONVECTION OF HEAT IN Liguips. Roy A. Nelson. Vol. 23, 
94-103, January, 1924.—On p. 101 the statement that A. H. 
Davis in his article in Phil. Mag. 44, 920, 1922, had assumed the heat 
loss per unit length to be proportional to @, is incorrect. In fact his 
formula No. 18 indicates that m lies between the limits given by the 
author. The mistake, which was pointed out by Mr. Davis, does not 
affect the conclusions of the paper but did an injustice to Mr. Davis 
which is sincerely regretted. 


A graphical method for the utilization of rotation spectra in crystal 
structure determination. Maurice L. Huggins. Vol. 23, 663, ab- 
stract 4, May, 1924,—The first equation should read 


i sin AV (h? /a®+k?/b*) —{n?\*(h?/a?+k?/b°+P/c*)? | 





1—4n?d?(h?/a?+k2/b?+2/c2) 








